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. Abstract 

D ' We study the dynamics of the restricted planar three-body problem near a mean motion 

, resonance, i.e. a resonance involving the Keplerian periods of the two lighter bodies revolving 

' around the most massive one. This problem is often used to model Sun- Jupiter-asteroid 

T-H , systems. For the primaries (Sun and Jupiter), we pick a realistic mass ratio fj, = 10~^ 

' and a small eccentricity eo > 0. The main result is a construction of a variety of diffusing 

, orbits which show a drastic change of the osculating eccentricity of the asteroid, while the 

' osculating semi major axis is kept almost constant. The proof relics on the careful analysis of 

. , the circular problem, which has a hyperbolic structure, but for which diffusion is prevented 

1^ ' by KAM tori. In the proof we verify certain non-degeneracy conditions numerically. 

, Based on the work of Treschev, it is natural to conjecture that the time of diffusion for 

^ ' this problem is at least ^ — ln(/xeo)/(/^'^/^eo). We expect our instability mechanism to apply 

I— I, to realistic values of eg and we give heuristic arguments in its favor. If so, the applicability 
of Nekhoroshev theory to the three-body problem as well as the long time stability become 

^ , questionable. 

■ It is well known that, in the Asteroid Belt, located between the orbits of Mars and 

, Jupiter, the distribution of asteroids has the so-called Kirkwood gaps exactly at mean motion 

OO ' resonances of low order. Our mechanism gives a possible explanation of their existence. To 

i relate the existence of Kirkwood gaps with Arnold diffusion, we also state a conjecture on 

On ' its existence for a typical e-pcrturbation of the product of the pendulum and the rotator. 

, Namely, we predict that a positive conditional measure of initial conditions concentrated in 

• the main resonance exhibits Arnold diffusion on time scales — Ine/e^. 
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1 Introduction and main result 



1.1 The problem of the stability of gravitating bodies 

The stability of the Solar System is a longstanding problem. Over the centuries, mathematicians 
and astronomers have spent an inordinate amount of energy proving stronger and stronger 
stability theorems for dynamical systems closely related to the Solar System, generally within 
the frame of the Newtonian A^-body problem: 

gv=^m,^^^^, (?iGR2, i = 0,l,...,N-l, (1) 

and its planetary subproblem, where mo (thought of as that of the Sun) is much larger than the 
other masses mj. 

A famous theorem of Lagrange entails that the observed variations in the motion of Jupiter 
and Saturn come from resonant terms of large amplitude and long period, but with zero average 
(see |Las06j and references therein, or |AKN88l Example 6.16]). Yet it is a mistake, which 
Laplace made, to infer the topological stability of the planetary system, since the theorem deals 
only with an approximation of the first order with respect to the masses, eccentricities and 
inclinations of the planets |Lap89 p. 296]. Another key result is Arnold's theorem, which proves 



the existence of a set of positive Lebesgue measure filled by invariant tori in planetary systems, 
provided that the masses of the planets are small |Arn63l Fej04| . However, in the phase space 



the gaps left by the invariant tori leave room for instability. 

It was a big surprise when the numerical computations of Sussman, Wisdom and Laskar 
showed that over the life span of the Sun, or even over a few million years, collisions and 
ejections of inner planets are probable (because of the exponential divergence of solutions, only 
a probabilistic result seems within the reach of numerical experiments); see for example |SW921 
ILas94j . or [LaslOj for a recent account. Our Solar System, as well as newly discovered extra-solar 
systems, are now widely believed to be unstable, and the general conjecture about the A^-body 
problem is quite the opposite of what it used to be: 

Conjecture 1.1 (Global instability of the A^-body problem). In restriction to any energy level 
of the N-body problem, the non-wandering set is nowhere dense. (One can reparameterize orbits 
so as to have a complete flow, despite collisions.) 

According to Herman |Her98j . this is the oldest open problem in dynamical systems (see 
also |Kol57j ). This conjecture would imply that bounded orbits form a nowhere dense set and 
that no topological stability whatsoever holds, in a very strong sense. It is largely confirmed 
by numerical experiments. In our Solar System, Laskar for instance has shown that collisions 
between Mars and Venus could occur within a few billion years. The coexistence of a nowhere 
dense set of positive measure of bounded quasiperiodic motions with an open and dense set of 
initial conditions with unbounded orbits is a remarkable conjecture. 

Currently the above conjecture is largely out of reach. A more modest but still much chal- 
lenging goal is a local version of the conjecture: 

Conjecture 1.2 (Instability of the planetary problem). If the masses of the planets are small 
enough, the wandering set accumulates on the set of circular, coplanar, Keplerian motions. 

There have been some prior attempts to prove such a conjecture. For instance, Moeckel 
discovered an instability mechanism in a special configuration of the 5-body problem |Moe96| . 



3 



His proof of diffusion was limited by the so-called big gaps problem between hyperbolic invariant 
tori; this problem was later solved in this setting by Zheng |ZhelO| . A somewhat opposite 
strategy was developed by Bolotin and McKay, using the Poincare orbits of the second species 
to show the existence of symbolic dynamics in the three-body problem, hence of chaotic orbits, 
but considering far from integrable, non-planetary conditions; see for example jBolOGj . Also, 
Delshams, Gidea and Roldan have shown an instability mechanism in the spatial restricted 
three-body problem, but only locally around the equilibrium point Li (see jPGRllj ). 

In this paper we prove the existence of large instabilities in a realistic planetary system 
and describe the associated instability mechanism. We thus provide a step towards the proof 
of Conjecture 11.21 The instability mechanism shown in this paper is related to a generalized 
version of the Mather mechanism |Mat96l IBT991 IDdlLSOOl IGT081 IKal03[ lPif06] . Some parts of 
the proof rely on numerical computations, but our strategy allows us to keep these computations 
to the simplest and, hopefully, most convincing. 

More specifically, consider the planetary problem ([1]), with one planet mass (say, mi), larger 
than the others small masses: mo S> mi ^ m2, ...,m7v-i- The equations of motion of lighter 
objects {i = 2, — 1) can advantageously be written as 



mo- 



^'"^^ Wn, -n.ii3 + 2^ _^.ii3 - 



\qo-qif wqi-qiW .^^^ "m-qi 



Letting the masses mj tend to for j = 2, A^ — 1, we obtain a collection of (A^ — 2) independent 
restricted problems: 

qo-qi , qi- qi 
\\qo-qi\r \\qi-qi\\ 

where the massless bodies are influenced by, without themselves influencing the primaries of 
masses mo and mi. 

For A^ = 3, this model is often used to approximate the dynamics of Sun-Jupiter-asteroid or 
other Sun-planet-object problems and it is the simplest one conjectured to have a wide range of 
instabilities. 

1.2 Relevance in astronomy 
1.2.1 The Asteroid Belt 

One place in the Solar system where the dynamics is well described by the restricted three- 
body problem, is the Asteroid Belt. The Asteroid Belt is located between the orbits of Mars 
and Jupiter and consists of 1.7 million objects ranging from asteroids of 950 kilometers to dust 
particles. Since the mass of Jupiter is approximately 2960 masses of Mars, away from close 
encounters with Mars, one can neglect the influence of Mars on the asteroids and focus on 
the influence of Jupiter. We also omit interactions with the second biggest planet in the Solar 
System, namely Saturn, which actually is not so small. Indeed, its mass is about a third of the 
mass of Jupiter and its semi major axis is about 1.83 times the semi major axis of Jupiter. This 
implies that the strength of interaction with Saturn is around 10% of the strength of interaction 
with Jupiter. However, instabilities discussed in this paper are fairly robust and we believe that 
they are not destroyed by the interaction with Saturn (or other celestial bodies), which, to some 
degree averages out. 

With these assumptions one can model the motion of the objects in the Asteroid Belt by 
the restricted problem. Denote by = mi/(mo + mi) the mass ratio, where mo is the mass of 
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the Sun and mi is the mass of Jupiter. For /i = (namely, neglecting the influence of Jupiter) 
bounded orbits of the asteroids are ellipses. Up to orientation the ellipses are characterized by 
their semi major axis a and eccentricity e. 

The aforementioned theorem of Lagrange asserts that, for /i > small, the semi major axis 
a{t) of an asteroid satisfies \a{t)—a{0)\ < for all \t\ < For very small n the time of stability 
was greatly improved by Niederman jNie96j using Nekhoroshev theory (see the discussion in the 
next section). Nevertheless, if one looks at the asteroid distribution in terms of their semi major 
axis, one encounters several gaps, the so-called Kirkwood gaps. It is believed that the existence 
of these gaps is due to instability mechanisms. 

1.2.2 Kirkwood gaps and Wisdom's ejection mechanism 

The mean motion resonances appear when the ratio between the period of Jupiter and the 
period of the asteroid is rational. The Kirkwood gaps correspond to the ratios 3:1, 5:2, 7:3. 



Asteroid Main-Belt Distribution 

Kirkwood Gaps 




Semi-major Axis (AU) 



Figure 1: Kirkwood gaps 



In this section we present a heuristic explanation of why these gaps exist. 

It is conjectured and confirmed by numerical data |Wis82j . that eccentricities of the asteroids 
appropriately placed in the Kirkwood gaps change by a magnitude of order of one. Notice that 
in the real data eccentricities of most asteroids in the Asteroid Belt are between and 0.25 (see 



e.g. http : //en . wikipedia . org/ wiki/File : Mainbeltevsa . png) . 

As the eccentricity of the asteroid grows while its semi major axis is nearly constant, its 
perihelion gets closer and closer to the origin, namely at the distance a{t)[l — e{t)], where a(t) 
and e{t) are the semi major axis and eccentricity of the asteroid respectively (see Figure [H 
where the inner circle is the orbit of Mars). In particular, a close encounter with Mars becomes 
increasingly probable. Eventually at some point Mars and the asteroid come close to one another 
and the asteroid most probably gets ejected from the Asteroid Belt. 

A surprising fact is that the change of eccentricity of the asteroid is only possible due to 
the ellipticity of the motion of Jupiter., due to the following count of dimensions. For circular 
motions of Jupiter the problem reduces to two degrees of freedom (see Section [2]) and plausibly 
there are invariant two-dimensional tori separating the three dimensional energy surfaces (see 
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e.g. GDF^SOt ICCOTj ). If the eccentricity of Jupiter is positive, the system has two and a half 



degrees of freedom and then KAM tori do not prevent drastic changes in the eccentricity. 

Heuristically, the conclusion is that, if the eccentricity of the asteroid changes by a magni- 
tude of order one in the Sun-Jupiter-asteroid restricted problem, then the asteroid might come 
into zones where the restricted problem does not describe dynamics appropriately, due to the 
influence of Mars. 

The main result of the paper is that in a certain mean motion resonance there are unstable 
motions which lead to significant changes in the eccentricity. In this paper we only present 
results for one particular resonance, the resonance 1 : 7, because the proof relies on numerical 
computations. This resonance is different from the ones in the Kirkwood gaps. It was initially 
chosen for reasons which turned out to be incidental (however, see the speculative remark in 
Section 11.4.31 about the speed of diffusion). Thus we are confident that our mechanism of 
instability applies to other resonances as long as the orbits of the unperturbed problem stay 
away from collisions. Thus, the instability mechanism showed in this paper gives an insight 
into the existence of the Kirkwood gaps. Another instability mechanism, using the adiabatic 
invariant theory, can be seen in jNS04j and applies to different ranges of mass ratio and Jupiter's 
eccentricity. 



1.2.3 Capture in resonance of other objects 

There are many known light objects in the Solar System which display a mean motion resonance 
of low order with Jupiter or some other planet. Some of them are: Trojan satellites (which 
librate around one of the two Lagrangian points of a planet, hence in 1 : 1 mean motion 
resonance with the planet), Uranus (which is close to the 7 : 1 mean motion resonance with 
Jupiter, thus giving an example of an "outer" restricted problem, close in phase space to the 
solutions we are studying), or the Kuiper Belt beyond Neptune (whose objects, behaving in the 
exact opposite manner to the thoses of the Asteroid Belt, seem to concentrate close to mean 
motion resonances; in particular, the Keplerian ellipse of the dwarf planet Pluto notoriously 
meets the ellipse of Neptune). The current existence of these resonant objects, and thus their 
relative stability, seemingly contradicts the above mechanism. This calls at least for a short 
explanation, although there are many effects at work here. 

The main point is that in the eye of a resonance lies an elliptic island, where some kind 
of long term stability prevails. Besides, the geometry often prevents the ejection mechanism 
described in Section 11.2.21 to occur in a straightforward manner (or to lead to a spectacular 
ejection), because there is no such body as Mars to kick the asteroid through a close encounter. 
In many cases, the mean motion resonance itself precludes collisions with the main planet (e.g. 
for Trojan asteroids with respect to Jupiter, obviously, or for Pluto with respect to Neptune; for 
a discussion of this effect in the Asteroid Belt, see |Rob05j ). One should add that the complete 
picture certainly includes secular resonances, close encounters between asteroids themselves, as 
well as more complicated kinds of resonance involving sometimes more bodies (e.g. the second 
Kirkwood gap, where a four-body problem resonance seems to play a crucial role). We refer to 
|Mor021 IRobOSj for further astronomical details. 



1.3 Main theorem 

Let us consider the three-body problem and let us assume that the massless body moves in the 
same plane as the pair of primaries. Normalize the total mass to one and call the three bodies 
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the Sun (mass 1 — fi), Jupiter (mass < /i ^ 1) and the asteroid (zero mass). If the energy 
of the primaries is negative, their orbits describe two elhpses with the same eccentricity, say 
Co > 0. For the sake of later simphcity, we wih change notations here, as we will denote by qo{t) 
the normalized position of the primaries ("fictitious body"), so that the Sun and Jupiter have 
respective positions —iJ,qo{t) and (1 — iJ,)qo{t). The Hamiltonian of the asteroid is 

1 1^ 1 1 ^ 1 — fi fj/ 

K{q,p,t) = ^ - J^^J^^ - (4) 

where q,p G M^. Without loss of generality one can assume that qo{t) has semi major axis equal 
to 1 and period 27r. For cq > this system has two and a half degrees of freedom. 

When Co = 0, the primaries describe uniform circular motions aroung their center of mass 
(restricted planar circular three-body problem). Hence, in a frame rotating with the primaries, 
the system becomes autonomous and thus has only 2 degrees of freedom. Its energy in the 
rotating frame is a first integral, called the Jacobi constant. It is defined by 

J = _ 1-/^ 

2 h + moim \\q - {I - fi)qo{t)\\ 2 ■ ^ ^ 

The aforementioned KAM theory applies to both the circular and the elliptic problems 
|Arn631 ISM95j and asserts that if the mass of Jupiter is small enough, there is a set of initial 
conditions of positive Lebesgue measure leading to quasiperiodic motions, in the neighborhood 
of circular motions of the asteroid. 

If Jupiter has a circular motion, since the system has only 2 degrees of freedom, KAM 
invariant tori are 2-dimensional and separate the 3-dimensional energy surfaces. But in the 
elliptic problem, 3-dimensional KAM tori do not prevent orbits to wander on a 5-dimensional 
phase space. In this paper we prove the existence of a wide enough set of wandering orbits in 
the elliptic planar restricted three-body problem. 

Let us write the Hamiltonian as 

K{q,p,t) = Ko{q,p) + Ki{q,p,t,fi) 

with 



Koiq,p) 
Ki{q,p,t,fi) 



\\pf 1 

2 ~ Ikll 

1 1 — A* ^J' 



\q + fJ-qoit)\\ lk-(l-^)go 



One can see that Ki = 0{fi) uniformly away from collisions. Then, notice that there is a 
competition between the integrability Kq and the non-integrability Ki, without which there 
would be no wandering. In this present work we consider a realistic value of the mass ratio, 
H = 10~^. For brevity in what follows we abbreviate the restricted planar circular (resp. elliptic) 
three-body problem to the circular (resp. elliptic) problem. 
Here is the main result of this paper. 

Theorem 1. Let us consider the elliptic problem with mass ratio fj, = and the eccentricity of 
Jupiter eo > 0. Then, for cq small enough, there exist T > and a trajectory whose eccentricity 
e{t) satisfies that 

e(0) < 0.48 and e{T) > 0.67. 
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Figure 2: Transition from the instant ellipse of eccentricty e = 0.48 to the ellipse of eccentricty 
e = 0.67. The dashed line schematically shows the transition. Nevertheless, the actual diffusing 
orbit is very complicated and diffusion is very slow. 



There is a more precise result, which is more difficult to state at this stage (see Theorem [5]). 
Let us say that along the above trajectory, the semi major axis a{t) remains almost constant, 
and, more precisely, 



(t) -7^/3 < 14/i for t G [0,r], 



and that 



|j(r)- j(o)| >o.i, 

where J is the Jacobi constant defined in ([5]). Now we discuss the main result from different 
perspectives. 



1.4 Refinements and comments 

1.4.1 Smallness of the eccentricity of Jupiter 

When Jupiter describes a circular motion, the Jacobi constant is an integral of motion and then 
KAM theory prevents global instabilities. We consider the eccentricity cq as a small parameter 
so that we can compare the dynamics of the elliptic problem with the dynamics of the circular 
one. 

It turns out that the difference between the elliptic and circular Hamiltonians is 0(/ieo). 
The analysis of the difference, performed later in the paper (see Section [4. 2p . shows that using 
averaging one can reduce the difference to ©(/iCq) (or even smaller). This makes us believe 
that eo should not be infinitesimally small for our mechanism to work. Even the actual value 
eo ~ 0.048 is not out of question. However, having realistic eo, to a large extend, becomes a 
matter of numerical experiment, not of mathematical proof -it is the limit and the interest of 
perturbation theory, to describe some dynamical behavior in terms of asymptotic models. See 
Appendix ID. 21 for more details. 
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1.4.2 On infinitesimally small /i's 



In Theorem [H we do not know what would happen asymptotically if we let — )• (our estimates 
worsen). Indeed, one of the crucial steps of the proof is to study the transversality of certain 
invariant manifolds (see Section ll.6p and this transversality becomes exponentially small with 
respect to ^ as — )• 0. On the other hand, Theorem [1] holds for realistic values of which is out 
of reach of many qualitative results of perturbation theory where parameters are conveniently 
assumed to be as small as needed. See Appendix ID] for more details. 



1.4.3 Speed of diffusion 

In Appendix |D] we discuss the relation of this problem with a priori unstable systems and the 
Mather accelerating problem. This leads us to conjecture that, for orbits constructed in this 
paper, the diffusion time T can be chosen 

T--'^. (6) 

De la Llave |dlL04j , Gelfreich-Turaev |GT08j , and Piftankin [PifOGj , using Treschev's techniques 
of separatrix maps (see e.g. |PT07j ). proved linear diffusion for the Mather acceleration problem. 
This might lead to possibility of linear diffusion after a smart choice of diffusing orbits, i.e. 
T ~ — ln/i(^'^/^eo)~^ (see Appendix iDl for more details) 0. 

We emphasize that an analytic proof might well require restrictive conditions between fi and 
eo. However, for realistic values of fi and eo or smaller, i.e. < n < 10~^ and < eo < 0.048, 
we expect that the speed of our mechanism of diffusion should also obey the above heuristic 
formula. 

On the other hand, the above asymptotics for the speed of diffusion probably does not hold 
in the neighborhood of circular motions of the masless body, which might be much more stable 
than more eccentric motions. This fact could give an explanation of the fact that Uranus, whose 
eccentricity of 0.04 is significantly smaller than most asteriods from the asteriod belt, and which 
is roughly in 7 : 1-resonance with Jupiter (its period is 7.11 times larger than that of Jupiter) 
has not been expelled yet (see also Section ri.2.3p . However, a deeper analysis would require to 
compare the distances of the various celestial bodies to the mean motion resonance, as well as 
the splitting of their invariant manifolds. 



1.4.4 On Nekhoroshev's stability 

Consider an analytic nearly integrable system of the form H^(9,I) = Hq(I) + eH\{6,I) with 
6' G and / in the unit ball B^. Suppose Hq is convex (or even weaker so-called steep). El 
Then a famous result of Nekhoroshev (see e.g. |Nie96j for the history and precise references and 
[XuelOj for the estimate on the involved constant c) states that for some c > independent of 
£ we have 

\I{t) - im < for |t| < exp(ce-V2n)^ 

Niederman |Nie96j applied Nekhoroshev theory to the planetary A^-body problem. He showed 
that the semi major axis obey the above estimate for exponentially long time exp(ce~^/^"') with e 

^This does not seem crucial, since the real value eo is not smaller than ^ 

■^Recall that Hq is called steep if for any afhne subspace L of R" the restriction Hq\l has only isolated critical 
points 
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being the smallness of the planetary masses. However, the constant c along with other constants 
involved in the proof are not so good. In particular, in order to have stability time comparable 
to the age of the Solar system e needs to be 3 • 10~^^. Moreover, having stability of semi major 
axis does not provide stability of the eccentricity, which we claim can have substantial deviations 
in polynomially long time. 

With e ~ /U, there was a hope that this result could be applied to long stability of e.g. the 
Sun-Jupiter-Saturn system (see e.g. |GG85j ). However, ([6]) indicates absence of even 0{e~'^)- 
stability. Indeed, the unperturbed Hamiltonian of the three body problem is neither convex, nor 
steep. This turns out to be not only a technical problem but a true obstruction to exponentially 
long stability, since Nekhoroshev's theory does not to apply to this kind of systems. See Appendix 
[D]for more details. 

1.5 Mechanism of instability 

The result obtained in Theorem [1] gives an example of large instability for this mechanical 
system. It can be interpreted as an example of Arnold diffusion (see [Arn64j ) . Nevertheless, 
Arnold diffusion usually refers to nearly integrable systems whereas Hamiltonian (jj]) cannot be 
considered as close to integrable since fi = is fixed. The mechanism of diffusion used in 
this paper is somewhat similar to the so-called Mather accelerating problem ( |Mat961 IBT991 
IDdlLSOOl IGT081 IKal03l [PifD6] 1. This analogy will be specified in Section E31 

Arguably, the main source of the existence of instabilities are resonances. One of the most 
natural resonances in the elliptic problem (even a three-body problem) are the mean motion 
orbital resonance]^. Along such a resonance, Jupiter and the asteroid will regularly be in the 
same relative position. Over a long time interval, Jupiter's influence could thus a priori pile up 
and, despite its small amplitude due to the small mass of Jupiter, could modify the eccentricity 
of the asteroid, instead of averaging out. According to Kepler's third Law, these resonances take 
place when a^/^ is close to a rational, where a is the semi major axis of instant ellipse of the 
asteroid. In our case we consider a^/^ close to 7. Nevertheless, one should expect that the same 
mechanism takes place for a large number of mean motion orbital resonances. 

The semi major axis a and the eccentricity e describe completely an instant ellipse of the 
asteroid (up to orientation). Therefore, geometrically Theorem [T] says that the asteroid evolves 
from a Keplerian ellipse of eccentricity e = 0.48 to one of eccentricity e = 0.67, without changing 
much its semi major axis (see Figure[2]). In Figure[3]we consider the plane (a, e), which describes 
the ellipse of the asteroid. Then diffusing orbits given by Theorem [1] correspond to a nearly 
horizontal line. 

A qualitative description of such a diffusing orbit is outlined at the end of section [5l 

1.6 Sketch of the proof 

The diffusing orbit of the elliptic problem we are looking for lies in a neighborhood of a (3- 
dimensional) normally hyperbolic invariant cylinder A and its local invariant manifolds, which 
exist near our mean motion resonance. The vertical component of the cylinder can be param- 
eterized by the eccentricity of the asteroid and the horizontal ones by its mean longitude and 
time. 

''The mean motions are the frequencies of the Keplerian revolution of Jupiter and the asteroid around the Sun: 
in our case the asteroid makes one full revolution while Jupiter makes seven revolutions. 
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Figure 3: In this graphic we show the diffusion path that we study in the (a, e) plane. The 
horizontal lines represent the resonances along which we drift. The thicker line is the diffusion 
path whose existence we are able to prove in this paper. 

If the stable and unstable invariant manifolds of A intersect transversally, the elliptic problem 
induces two different dynamics on the cylinder (see Sections 14.41 and 14. 5|) : the inner and the 
outer ones. The inner dynamics is simply the restriction of the Newtonian flow to A. The 
outer dynamics is obtained by a limiting process: it can be observed asymptotically by starting 
very close to the cylinder and to its unstable manifold, traveling all the way up to a homoclinic 
intersection, and coming back close to the cylinder and along its stable manifold. Since the 
system has different homoclinic orbits to the cylinder, one can define several different outer 
dynamics. In our diffusing mechanism we use two different outer maps. The reason is that each 
of the outer maps are not defined in the whole cylinder and then we need the two of them to 
achieve diffusion (see Section [3]). 

The proof consists in the following four steps: 

1. Prove existence of the normally hyperbolic invariant cylinder A. 

2. Establish transversality of the stable and unstable invariant manifolds of this cylinder. 

3. Compare the inner and outer dynamics on A and, in particular, check that they do not 
have common invariant circles. 

4. Construct diffusing orbits by shadowing a carefully chosen composition of the outer and 
inner maps. 

This program faces difficulties at every step. 

1.6.1 Existence of a normally hyperbolic invariant cylinder A 

The first difficulty comes from the proper degeneracy of the Newtonian potential: at the limit ^ = 
(no Jupiter), the asteroid has a one- frequency, Keplerian motion, whereas symplectic geometry 
would allow for a three-frequency motion (as with any potential other than the Newtonian 
potential 1/r and the elastic potential r^). Due to this degeneracy, switching to /i > (even 
with eo = 0) is a singular perturbation. 
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1.6.2 Transversality of invariant manifolds 

The second step, establishing the transversahty of the invariant manifolds of A, even for cq = 0, is 
a delicate problem. Asymptotically when ^ ^ 1, the difference (splitting) between the invariant 
manifolds becomes exponentially small with respect to fi, that is of order exp(— c/y^) for some 
constant c > 0. Despite inordinate efforts of specialists, all the known techniques fail to estimate 
this splitting, because the relevant Poincare-Melnikov integral is not algebraic. Note that this 
step simplifies dramatically when people study generic systems. 

At the expense of creating other difficulties, setting fi = 10~^ avoids this splitting problem, 
since for this value of the parameter one can see that the splitting of separatrices is not extremely 
small and therefore, can be detected by means of a computer with convincing accuracy. Besides, 
= 10~^ is a realistic value of the mass ratio for the Sun-Jupiter model. Since the splitting of 
the separatrices varies smoothly with respect to the eccentricity eg of the primaries, it suffices 
to estimate the splitting numerically for cq = 0, i.e. in the circular problem. This is a key point 
for the numerical computation, which thus remains relatively simple. On the other hand, in the 
following two steps it will be crucial to have eo > 0, without which the KAM tori would separate 
energy levels. 

Finally, recall that the cylinder A has two branches of both stable and unstable invariant 
manifolds. In certain regions the intersections between one of the branches of the stable and 
unstable invariant manifolds is tangential, which does not allow to define the outer map. Nev- 
ertheless, then one can check that the other two branches intersect transversally so that we 
can define a different outer map. Thus, we will combine the two outer maps depending which 
branches of the invariant manifolds intersect transversally. 



1.6.3 Asymptotic formulas for the outer and inner maps 

Now we turn to the third step. Using classical perturbation theory and the specific properties of 
the underlying system one can reduce the inner and (the two different) outer dynamics to three 
two-dimensional symplectic smooth maps of the form 

■ V * y V i + /uro(/,/i) + 0(/ieo) 

and 

^out,* ( / + eo(i?*'+(/,^i)e»* + i?*'"(/,^)e-'*)+0(/ie2) 

where (/, t) are conjugate variables which parameterize a certain connected component of the 
3-dimensional normally hyperbolic invariant cylinder A intersected with a certain transversal 
Poincare section and , To, B*'^ ,uj* are smooth functions. The superindexes f and b stand for 
the forward and backward heteroclinic orbits that are used to define the outer maps. The choice 
of this notation will be clear later on in Section [3l Note that these maps are real and therefore 
A~ and i?*'^ are complex conjugate to A'^ and B*'~^ respectively. 




* = f,b, (8) 



1.6.4 Non-degeneracy implies the existence of diffusing orbits 

As we will see in Section \5\ the existence of diffusing orbits can be established provided the 
smooth functions 
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do not vanish for all / G [/_,/+] for which the corresponding outer map is defined. Since and 
B*'^ are complex conjugate, we do not write the complex conjugate /C*'~(/, /i). Numerically, 
one can check that 1C*'~^(I, /i) 7^ in their domain of definition. It turns out that fC*'~^(I, /x) 7^ 
implies absence of common invariant curves for the inner and outer maps. This reduces the 
proof of Theorem [1] to shadowing, made in step 4, and thus it leads to the existence of diffusing 
orbits. Moreover, it turns out that for this problem no large gaps appear. This fact is not so 
surprising taking into account that the elliptic problem has three time scales. 

Finally, let us point out that the complex functions IC*'~^{I,n) can be regarded as a 2- 
dimensional real- valued function depending smoothly on (/,//). If the dependence in /i is non- 
trivial, a complex valued function /C*'"^ (/,//) does not vanish at any point of their domain of 
definition except for a finite number of /x's. 

1.7 Nature of numerics 

In this section we outline which parts of the mechanism are based on numerics. 

• On each 3-dimensional energy surface the circular problem has a well-defined Poincare 
map Fj : Sj — )■ Sj of a 2-dimensional cylinder Sj for a range of J's. For each J in some 
interval [J_, J-|_] we establish the existence of a saddle periodic orbit Fj{pj) = pj. 

• We show that for all J E [J_, J+] we have two intersections of W^{pj) and W^{pj). 
Each intersection is transversal for almost all values of J, but it becomes tangent at an 
exceptional (discrete) set of J's. Nevertheless, we check that at least one of the two 
intersections is transversal for each J G [J_, J^,]. See figure [TCI 

• Each transversal intersection qj gives rise to a homoclinic orbit, denoted jj. For each 
J G [ J_ , J+] we compute several Melnikov integrals of certain quantities related to AHdi 
along 7 J and pj. Out of these integrals we compute the leading terms of the dynamics of 
the elliptic problem and verify a necessary condition for diffusion. 

As seen in the Appendixes ElCl the numerical values that we deal with are several orders of 
magnitude larger than the estimated error of our computations, and therefore these computations 
are reliable. Moreover, all the computations that we handle are standard and low-dimensional. 

2 Setting of and notations 

The model of the Sun, Jupiter and a massless asteroid in Cartesian coordinates is given by the 
Hamiltonian ([4]). First, let us consider the case /i = 0, that is, we consider Jupiter with zero 
mass. In that case, Jupiter and the asteroid do not make influence on each other and therefore 
the system is reduced to two uncoupled 2-body problems, the Sun-Jupiter and the Sun-asteroid, 
which are integrable. 

We want to study the existence of instability in one particular resonance of this system, 
which appears when the period of the asteroid is seven times the period of Jupiter. One can 
consider the so-called Delaunay variables, which we denote by (£, L, g, G), which are angle-action 
coordinates of the Sun-asteroid system. The variable £ is the mean anomaly, L is the square 
of the semi major axis, g is the argument of the perihelion and G is the angular momentum. 
These variables can be obtained from the Cartesian coordinates as follows (see |AKN88) for more 
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details and background, or |FejlO[ Appendix] for a straightforward definition). First define polar 
coordinates for the position: 

q = (r cos(/>, rsincj)). 



Then, the actions of the Delaunay coordinates are defined by 

(10) 



1 1 



2L2 2 



\Q\ 



= -J-^ (11) 

(recall that /i = for these definitions). Using these actions, the eccentricity of the asteroid can 
be expressed as 

1 - JJ- (12) 

To define the angles, let v and g be the true anomaly and the argument of the perihelion, so 
that 

^ = v + g. (13) 
Then, from v one can obtain the eccentric anomaly u using 



V 1 + e u 
tan- = y^tan-. (14) 

From the eccentric anomaly, the mean anomaly is given by Kepler's equation 

u — esinti = i. (15) 

In Delaunay coordinates, the Hamiltonian ^ can be split into the Keplerian part —1/2L'^, 
the circular part of the perturbing function fiAHdrc and the remainder which vanishes when 
eo = 0: 

if(L, e,G,g-t,t) = -^ + /iAifcirc(^, i,G,g- t, fi) + fieoAH,n{L, e,G,g- t, t, fi, eo). (16) 

For eo = 0, the circular problem only depends on g — t. To simplify the comparison with the 
circular problem, we consider rotating Delaunay coordinates, in which AHdrc is autonomous. 
Define the new angle g = g — t (the argument of the pericenter, measured in the rotating frame) 
and a new variable / conjugate to time t. Then, we have 

H{L, l,G,gJ,t) = -^-G + ^iAH^i,c{L, £, G, g, fi) + fieoAH^n{L, i, G, g, t, fi,eo)+ I. (17) 

In these new variables, the difference of number of degrees of freedom of the elliptic and circular 
problems becomes more apparent. When eo = the system is autonomous and then / is constant, 
which corresponds to the conservation of the Jacobi constant ([5]). Therefore, the circular problem 
reduces to 2 degrees of freedom. Moreover, it will later be crucial to see the circular problem 
as an approximation of the elliptic one, in order to reduce the Herculean (and more doubtful) 
numerical computations of a direct approach to the corresponding lower dimensional, and thus 
simpler, computations of the circular problem. 
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Recall that we consider the 1 : 7 mean motion orbit resonance between Jupiter and the 
asteroid. That is, the period of the asteroid being approximately seven times the period of 
Jupiter. In rotating Delaunay variables, this corresponds to 

and g 1. (18) 

A nearby resonance is 

£ ~ - and t ~ 1, 
7 

but we will stick to the previous one. 

The resonance takes place when L ~ 7^^'^. We will study the dynamics in a large neigh- 
borhood of this resonance and we will see that one can drift along it. Namely, we will find 
trajectories that keep L close to 7^^^ while the G-component changes noticeably. Using ()12p . 
one can see that e also changes by an order of one. In this setting. Theorem [1] can be rephrased 
as follows. 

Theorem 2. There exist Cq > such that for < cq < Eq, there exist T > and an orbit of 
the Hamiltonian System with Hamiltonian (fT7|) which satisfies 

G(0) < Go and G{T) > d 

whereas 



L{t) - 7^/3 



< 7/i, 



By definition the Hamiltonian (jl7p is autonomous and thus preserved. Therefore, we will 
restrict ourselves to a level of energy which, without loss of generality, can be taken as = 0. 
Therefore, since |/ — G| = 0{fi), for orbits satisfying \L{t) — 7^/^| < 7fi, drift in G is equivalent 
to drift in /. 

The proof of this theorem is structured as follows. 

In Sectional we study the dynamics in the circular problem, that is eo = and the underlying 
Hamiltonian ()17p becomes 

H,i,,{L,£,G,g) = -^-G + fiAH,i,,{L,i,G,g,fi). (19) 

Theorem ^ says that for an interval of Jacobi energies [J_,J+] the circular problem has a 
smooth family of hyperbolic periodic orbits Aj, whose stable and unstable manifolds intersect 
transversally for each J G [J_, J+]. This theorem implies (Corollary [3Tj existence of a normally 
hyperbolic invariant cylinder. Later in the section (Subsections 13.21 and 13. 3p we calculate the 
aforementioned outer and inner maps for the circular problem (see ([7]) and ([8])). 

Then in Section H] we consider the elliptic case eo > as a perturbation of the circular case. 

Theorem 2] says that the family of periodic orbits give rise to a normally 

hyperbolic invariant cylinder A^q whose stable and unstable manifolds intersect transversally 
for each J G [J_ + 6, J+ — 6] with small 6 > 0. These objects give rise to the inner and outer 
maps for the elliptic problem. T/ieorern [5] provides expansions for the inner and outer maps (see 
formulas (|46p and (|49p respectively). 

Finally, in Section [5] in Theorem [6] we complete the proof of Theorem [2j This is done by 
comparing the inner and the two outer maps in Lemma 15.21 and constructing a transition chain 
of tori. It turns out that there are no large gaps, due to the specific structure of times scales 
and the Fourier series involved. This a priori contrasts with the typical situation of dynamics 
near a resonance (see e.g. |DdlLS06] ). 
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Notation 2.1. From now on, we will omit the dependence on fi (keeping in mind at various 
points the question of what would happen if we let fj, vary). Recall that we are taking a realistic 
value of fi = 10~^. 



3 The circular problem 
3.1 Invariant cylinder 

The circular problem is given by the Hamiltonian (|17p with eo = 0. Since it does not depend 
on t, I is an integral of motion. Moreover, since we are studying the dynamics in the energy 
surface H = 0, we have / = —HciTc{^,L,g,G). Therefore, the variable / equals the opposite of 
the Jacobi constant (0). For each level / = constant, one can study the dynamics close to the 
resonance 7i + g ^ 0. Since t is a cyclic variable, one can consider the two degree of freedom 
Hamiltonian of the circular problem for which the conservation of energy corresponds to the 
conservation of the Jacobi constant ([5]). Moreover, one can see that the circular problem is 
reversible with respect to the involution 

^{L,i,G,g,I,t) = {L,-e,G,-g,I,-t). (20) 

This symmetry will facilitate several numerical computations. 

Theorem 3. Consider the Hamiltonian (jl9p with fi = 10"'^. Then, in each energy level J G 
[J_, J+] = [—1.81,-1.56], there exists a hyperbolic periodic orbit Xj = {Lj{t),ij{t),Gj{t),gj{t)) 
of period Tj which satisfies 

\Tj - Utt\ < 60fi, 



and is smooth with respect to J, and 



Lj{t) - 7^/3 



< 7n 



for all t G M. 

Each Xj has two branches of stable and unstable invariant manifolds W^'^{Xj) and W^'^{Xj), 
j = 1,2. Then, for each J G [J_, J+] either W''^{Xj) and W'^iXj) or W'''^{Xj) and W'^iXj) 
intersect transver sally. 

Proof. Based on convincing numerical data. See Appendix |X1 □ 

We will study the elliptic problem as a perturbation of the circular one. Therefore we 
do not reduce the dimension of the phase space while studying the inner and outer dynamics 
of the circular problem. Namely, we consider the Extended Circular Problem given by the 
Hamiltonian ()17p with eo = 0. In other words, we keep the conjugate variables {I,t) even if t 
is a cyclic variable. Consider the energy level H = 0. In this setting the conservation of the 
Jacobi constant corresponds to the conservation of /. Therefore, the periodic orbits obtained 
in Theorem [3] become invariant two-dimensional tori which belong to hyperplanes / = constant 
for any 

/G [/_,/+] = [-J+,-J_] = [1.56,1.81]. (21) 

Moreover, the union of these 2-dimensional invariant tori form a normally hyperbolic invariant 
3-dimensional manifold. 
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Corollary 3.1. The Hamiltonian ()17p with fi = 10~^ and eg = has an analytic normally 
hyperbolic invariant 3 -dimensional manifold Aq, which is foliated by two-dimensional invariant 
tori. 

Moreover, Aq has two branches of stable and unstable invariant manifolds, which we call 
W^'^{Aq) and W'^'^{Aq), j = 1,2. Then, in the invariant planes I = constant, for each I G 
[/_,/+] either W'^'^{Aq) and W^'^{Aq) orVF'^'^(Ao) and VF"'^(Ao) intersect transversally. 

It follows from the implicit function theorem applied with the energy as a parameter (the 
nondegeneracy of the hyperbolic periodic orbits being observed numerically) that the cylinder 
Aq is analytic. 

We define a global Poincare section and deal with maps to reduce the dimension by one. 
There are two natural choices: {t = 0} and {g = 0}, since both variables satisfy t / and g ^ 0. 
We choose the section {g = 0} and call 

VQ:{g = 0} ^{g = 0} (22) 

this Poincare map. Since we are studying the resonance (jlSp . the intersection of the cylinder Aq 
with the section = 0} is formed by seven cylinders (see Figure H]). We denote them by Aq, 
j = 0, . . . , 6. Namely, 

_ 6 _ 

Aon{5 = 0}=Ao = U^o- (23) 

j=0 

As a whole U^^qAq is a normally hyperbolic invariant manifold for the Poincare map Vq. One 
can also consider the Poincare map Vq, namely iterate seven times Vq. Then, for this map, 
each Aq is a normally hyperbolic invariant manifold (of course, their union is also a normally 
hyperbolic invariant manifold for Vq). We work with the cylinders Aq since they have the 
advantage of having a natural system of coordinates. This system of coordinates will be used 
later on to study the inner and outer dynamics on them. In particular, we will work with Aq 
and A^. The reason is that in each invariant plane / = constant they are connected by at least 
one heteroclinic connection (of Vq) which is symmetric with respect to the involution ()20p . We 
call it a forward heteroclinic orbit if it is asymptotic to Aq in the past and Aq in the future and 
a backward heteroclinic orbit if it is asymptotic to Aq in the past and to Aq in the future. We 
denote by C [/_,/+], where f stands for forward, the subset of where W^{Aq) and 

M^*(Aq) intersect transversally and by C [/_,/+], where b stands for backward, the subset 
of [/_,/+] where I^'^(Aq) and W^{Aq) intersect transversally. By Corollarv 13.11 we have that 
pf UP^ = [/_,/+]. 

Corollary 3.2. The Poincare map Vq defined in (|22p . which is induced by the Hamiltonian 
(|17p with fj, = 10~^ and cq = 0, has seven analytic normally hyperbolic invariant manifolds 
Aq, j = 0, . . . ,6. They are foliated by one- dimensional invariant curves. Moreover, there exist 
analytic functions : [/_,/+] x T ^ (M x T)^ , 

Gl{I,t) = [Gi{I),o,I,t) = {g^^\l),G^^\l),G=f{I),Q,I,t) , (24) 

that parameterize Aq, namely, 

A^Q = {g^(I,t):(I,t)G [/„,/+] xt}. 



17 



Figure 4: The periodic orbit obtained at each energy level intersects seven times the Poincare 
section {g = 0}, as it is shown schematically in this picture. Then, when one considers the 
Poincare map Vq, the normally hyperbolic invariant manifold Aq has seven connected compo- 
nents Aq,. . . , Aq- 

The associated invariant manifolds VF"(Aq) and W^{Aq) intersect transversally provided I G ; 
and VF*(Aq) and W^{Aq) intersect transversally provided I G 2?^. Moreover, one of the points 
of these intersections belongs to the symmetry axis of (j20p . Let us denote by Tq, * = f,b, these 
intersections. Then, there exist analytic functions 

C(5 xM^ (MxT)^ {I,t) ^ C^{I,t), * = f,b 

which parameterize them: 

PS = [Q{I,t) = (C*'^(/),C*''(/),C*'^(/),0,/,t) : {I,t) G X t} , * = f,b. 

The zero which appears in the parameterizations G and C is just the ^(-coordinate. We keep 
it although we are in the Poincare section because later on we will use these parameterizations 
in the full phase space. 

Again, it follows from the implicit function theorem that M^*(Aq) and VK"(Ag) are analytic 
(taking the distance from the cylinder Aj] or Aq as a small parameter, as in |Mey75| with the 
cylinder as factor variable). 

Corollarv 13.11 gives global coordinates (/, t) for each cylinder Aq. These coordinates are 
symplectic with respect to the canonical symplectic form 

no = dI A dt, (25) 

Indeed, one has to consider the pullback of the canonical form dL Adi + dG Adg + dl Adt to the 
cylinders A;j. By Corollarv 13. 21 we have that in the cylinders: g = ^, t = Q^'^{I) and L = Gq^{I)- 
Then, it is easy to see that the pullback of dL A d£ + dG Adg + dL Adt is just VLq. 

We consider the inner and the two outer maps in one of these cylinders. We choose Aq. 
As we have explained before, the reason is that the heteroclinic connections with the following 
cylinder Aq intersect the symmetry axis of the involution (|2U|) and therefore they are easier to 
be studied numerically (see Figure [TT]) . Since I is conserved by the inner and outer maps, these 
maps are integrable and the variables (/, t) are the action-angle variables. In this way, it will be 
easier to understand the influence of ellipticity. 
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3.2 The inner map 



We first define the inner map. One could take the Poincare map Vq restricted to the whole 
normally hyperbolic invariant manifold ()23p . Nevertheless, in order to study the diff'usion mech- 
anism it is more convenient to consider just one of the cylinders that form (|23p . for instance 
Aq. To this end, we can define the inner map J-"™ : Ag — >• Aq as the analytic Poincare map 
Vq restricted to the symplectic invariant submanifold Aq. We can express J-"™ using the global 
coordinates (/, t) of Aq. 

Since / is an integral of motion, this inner map is of the form 

where the function To is independent of t due to the fact that the inner map preserves the 
differential form (|25p . which does not depend on t, and that / is a first integral. In fact, 
147r + fiTo{I) is the period of the periodic orbit obtained in Theorem [3] on the corresponding 
energy surface. It can be seen numerically that this map is twist. 

Lemma 3.3. The analytic symplectic inner map J-^ defined in (j26p is twist, that is 

diTo{I)^0 for I € [1^,1+]. 

Moreover, the function To{I) satisfies 

< fiTo{I) < TT. (27) 
Proof. Based on convincing numerical data. See Appendix |Al □ 



In Section [331 the function To{I) will be written by means of an integral (see ([39]) ). 
The information contained in this lemma will be crucial in Section [5] to prove the existence 
of a transition chain of invariant tori. 



3.3 Outer map 

Sometimes the outer map is called scattering map (see for instance |Tre04[ [DdlLSOS] ). We first 
recall the construction in a general perturbative setting, which we apply next to the circular 
problem, and in section 14.11 to the elliptic problem. 

Let Vq be a map of a compact manifold M. Let Aq C M be a normally hyperbolic invariant 
manifold ofVo, whose inner map PoUo ^^ro Lyapunov exponents: lim„_5.+oo In \\dVQ{z)v\\/n = 
for any z £ Aq and v G T^Aq (where || • || is some smooth Riemannian norm on M). Further 
assume that the stable and unstable invariant manifolds of Aq intersect transversally. 

Let V he a small perturbation of Vq. Since Aq is normally hyperbolic it persists under small 
perturbation of Vq. Let A C M be a normally hyperbolic invariant manifold of V. 

Then, the outer map associated to V and A (a particular case being V = Vq and A = Aq) is 
defined over some domain as follows. 

Definition 3.4. Assume that 7 C H WJ^ is a homoclinic orbit and that the intersection of 
and Wj( is transversal along 7, that is 

T,Wi + T,Wl = T,M and T,Wlr\T,Wl = T,^ forzG^. 

Then, we say that 5(x_) =2;+, if there exists a point z G 7 such that for some C > we have 

dist(P"(z),P"(x±)) < CA-I"I for all n G (28) 
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That condition ()28p indeed defines a map x- i— )• locally uniquely, is justified in |DdlLS08] . 
Remark 3.5. Since A is normally hyperbolic, for each point x £ A there are strong stable and 
unstable manifolds W^^{x) and W'^'^{x). Then = x+ holds only if VF*'"(x_)nVF*^(3;+) / 

and the intersection occurs on 7. 

In the case when the Lyapunov exponents of the inner dynamics VIa are positive, for the 
points x_ and x^ to be still uniquely defined given 2: G 7, A must exceed the maximal Lyapunov 
exponent i.e., the convergence towards A must dominate the motion inside of A. Otherwise, one 
cannot distinguish if the orbit of z is (backward- or forward-) asymptotic to a point of A or to 
the stable manifold of this point. 

Remark 3.6. If the Lyapunov exponents of the inner map V\a are zero (thus in particular for 
the unperturbed map Vq), the outer map S is C°°. 

On the other hand, if the Lyapunov exponents of the inner map are small (thus in particular 
for a map V close enough to Vq), the outer map is C^, and k tends to infinity as the Lyapunov 
exponents tend to 0. 

Strictly speaking, there is hardly any published regularity theorem from which these asser- 
tions follow directly. In order to prove them, one can first localize in the neighborhood of a 
small continuous set of hyperbolic periodic orbits of Vq, modify V outside this neighborhood 
in order to embed the periodic orbits into a compact invariant normally hyperbolic cylinder, 
and characterize for the stable and unstable manifolds of the modified system in terms of an 
equation of class C^, the perturbative parameter being the distance from the invariant cylinder. 
Such arguments belong to well understood hyperbolic theory, and we will omit further details, 
refering to the techniques developped in |Fen721 ICha04] or [BKZIH Appendix B] for a closer 
context. 

We now intend to apply a variant of this definition, in the unperturbed case, to the dynamics 
of the circular problem. One could consider the normally hyperbolic invariant manifold uJ^^Aq. 
Nevertheless, since it is not connected, it is more convenient to proceed as we have done for 
the inner map in the previous section. Namely, we look for an outer map which sends Aq to 
itself. Now one has to be more careful since one cannot consider the outer maps induced by the 
Poincare map Vq. Indeed, for Vq, the cylinder Aq is a normally hyperbolic invariant manifold 
but the homoclinic points obtained in Theorem 13.11 now correspond to heteroclinic connections 
between Ag and Aq and between Aq and Ag. To overcome this problem we will compose the 
heteroclinic outer maps 5*, * = f, b with the Poincare map Vq as many times as necessary so 
that the composition sends Aq to itself. 

Therefore, the smooth outer maps J^"*'^ that we consider connect Ag to itself and are 
defined as 



(29) 

Q - O u /-Q . JIQ r ilg. 



where is the outer map which connects Aq and Aq through VF"(Ag) n VF'^(Ag) and is the 
outer map which connects Aq and Ag through 14^" (Aq) n 14^'^ (Ag). Note that here we are abusing 
notation since the forward and backwards outer maps are only defined provided I G 2?^ and 
I G T)^ respectively and not in the whole cylinder Aq. 

The outer map is always exact symplectic (see [DdlLSOS] ). So, in the circular problem, since 
/ is preserved, the outer maps must be of the form 



^ M,j-(, + ^.(„).-=f.'=- (30) 
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Outer maps can be defined with either discrete or continuous time. Since the Poincare-Melnikov 
theory is considerably simpler for flows than for maps, we compute J^"*'* using continuous time. 
Moreover, in Section 14.51 we will use also flows to study the outer map of the elliptic problem 
as a perturbation of ([HH]) . The outer map given by the Hamiltonian ([T7D with cq = does 
not preserve the section {g = 0} but the inner map does. In order to fix this problem, we 
reparameterize the flow so that the inner and outer map preserve this section. 

This reparameterization corresponds to identifying the variable g with time and is given by, 

d „ BlH d diH 



ds —1 + fldc^Hcirc ds —1 + fidc^Hcirc 

— Q= I —G= (31) 

ds ds -1 + fidcAHcirc 

±t= ^ ^1= 

ds — 1 + fidcAHcirc ds 

where H is Hamiltonian (jl7p with cq = 0. Notice that this reparameterization implies the 
change of direction of time. However, the geometric objects stay the same. In particular, the 
new flow still possesses the normally invariant cylinder obtained in Corollarv 13 . 1 1 and its invariant 
manifolds. 

We will refer to this system as a reduced circular problem. We call it reduced because we 
identify g with the time s. Moreover, the t-component (and in fact, all the others) only depends 
on the other coordinates. Denote by ^q'"'^ the flow associated to the {L,i,G,g) components of 
equation (j31|) (which are independent of t and /). Componentwise it can be written as 

^^''{s,iLJ,G,g)} = (^<^(^{s,{LJ,G,g)},^'o{s,{L,£,G,g)},^^{s,{L,i,G,g)},g + s). (32) 

Then, the outer map can be computed as follows. 
Let 

J* (a) = cI>r{cT, (C*'^/),C*'^(/),C*'^(/),0)}, * = f,b 

Ai(a) = <i>r{a,igi'\i),gl,'\i),gfii),o)} 



be trajectories of the circular problem. The flrst ones have the initial conditions at the homo- 
clinic points obtained in Theorem [3] with action / since Cq ^ are the parameterizations of the 
intersections of the invariant manifolds of Aq and Aq , given in Corollary 13.21 The second ones 
have the initial condition in the flxed points of the Poincare map Vq, which are parameterized 
by Qq and given in Corollary 13.21 

Lemma 3.7. The functions uj^'^{I) involved in the definition of the outer maps in (|30|) can be 
defined as 

a;*(/) =<,,(/) +<(/), 



where 



with 



col^dl) = - ^1(1) (34) 



147V7r 



io^{I) = Jim / , , ,,r, . ^ — L.*r. da + N%{I) , * = f,b 



N^-oo \Jq -1 + l2{dG AHcirc) o 7| (o") 



(35) 
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and 



where To{I) is the function in (j26p . 

Note that the minus sign that appears in the hmits of integration of the integrals involved in 
the definition of functions is due to the fact that the reparameterized flow ()3ip reverses 

time. Using that the circular problem is symmetric with respect to ()20p and that the homoclinic 
points Cq and Cq belong to the symmetry axis, one can easily see that = —uj"^, * = f, b. 

The geometric interpretation of U!^'^{I) is that the t-shift occurs since the homoclinic orbits 
approach different points of the same invariant curve in the future and in the past. This shift 
is equivalent to the shift in t that appears in the Mather Problem |Mat96j . See, for instance, 
formula (2.1) in Theorem 2.1 of [DdlLSOO] and the constants a and b used in formula (1.4) of 
|BT99j . 

Proof. We compute uj^{I)- The other case can be done analogously. Since the t-component of 
the reduced circular system pip does not depend on t, its behavior is given by 



^l{s, (L, A G, g, t)}=t + $o{5, (L, £, G, g)} 

(37) 



= t + / ; — : r da 

Jo -1 + ^ScAifeirc {'^rW, {L,£, G,g)}) 

Note that, using this flow, the inner map on (|26p is just the (— 147r)-time map in the time s. 
Then, the original period of the periodic orbits obtained in Theorem [3l can be expressed using 
this new flow as 

Jq -1 + fl{dG AHcircjo Xjicr) 
This allows us to define the function To{I) in (|26p through integrals as 



-1 + nidcAH^irc) o Af (a) 

Consider now a point (Co'-^(I),Co'^(/), Co'^(I), 0, /, t) mW''{A.l)nW'(A^)n{g = 0}. Since the 
first four components are independent of t, (Cq '^(I), Cq^(I), Cq *^(/), 0, 1, t) is forward asymptotic 
(in the reparameterized time) to a point 



and backward asymptotic (in the reparameterized time) to a point 

fe'^(/), g'/{i), aQ''^(i), 0, /, t + . 
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Using ([37|) . the functions (^±{I) can be defined as 

1 1 



u}Ul)= lim [ I 

co^ (I) = lim [ I 



1 1 



+ fiidcAH^i,,) o 71 (a) -1 + /i(aGA/7circ) o A|(a) 



I da 

(40) 

da. 



Moreover, since the system is 147r-periodic in the time s due to the identification of s with g, it 
is more convenient to write these in integrals as 



LoUl) = lim / 

J (I) = lim / 




+ fiidoAH,;,,) o 7} (cj) -1 + fiidcAH,;,,) o A|(a) 
Then, taking (|38p into account, one obtains 



(/•MNtt 1 \ 

Jo -1 + fl{dGAHcirc) o 7/(cr) / 

from which the formulas for in ()35p follow. 

Finally we have to compute ujI^{I). This term corresponds to the contribution of Vq to 

the outer map in formula (|29p . Then, taking into account that t is defined modulo 27r is 

straightforward to obtain <^[n(/) in (|35p . □ 



4 The elliptic problem 



We can eventually study the elliptic problem. We will obtain perturbative expansions of the 
inner and outer maps. To this end, we use Poincare-Melnikov techniques applied to the reduced 
elliptic problem, which is given by 



ds* 



Ts9 



-1 + ndcAHarc + neodcAHeii 



ds 



L 



diH 



1 



-1 + jj-dcAHchc + i^eodcAHcw 



ds^ 



-1 + fidaAHavc + fieodaAHcu 



-1 + ^idcAHcnc + fJ-eodaAHeii 

-1 + ^idaAHcnc + fJ-eodaAHeii ' 



(41) 

This system is a perturbation of (|3ip . The study of the inner map can be done both dealing 
with this system or the system associated to the Hamiltonian (|16p . Nevertheless, to simplify the 
exposition we use only (|4ip for both the inner and outer maps. We also consider the Poincare 
map associated with this system and the section {g = 0}, 



Veo ■{9 = 0}^{g = 0}, 



(42) 



which is a perturbation of 

There are two main results in this section: 
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• Existence of a normally hyperbolic invariant manifold with transversal intersections of its 
invariant manifolds for the elliptic problem (Theorem [3D 

• Computation of the eo-expansions of the inner and outer maps associated to it (Theorem 

The above theorems are stated in the next section. Theorem |3] relies on Corollary[321 because 
we study the elliptic problem as a perturbation of the circular one. The proof of Theorem [5] 
consists of several steps. In Section [4.21 we analyze the eo-expansion of the elliptic Hamiltonian, 
and from it, in Section 14.31 we deduce properties of the eo-expansion of the flow associated to 
system (j41|) . In Section [4.41 we do perturbative analysis of the normally hyperbolic invariant 
cylinders Aig , which are the perturbation of the cylinders Aq obtained in Corollary 13.21 This 
analysis allow us to obtain perturbative in cq formulas for the inner map derived. Finally, in 
Section 14.51 we use the above expansions and compute the outer maps using Poincare-Melnikov 
techniques. These expansions allows us to give perturbative formulas of the inner and outer 
maps of the elliptic problem, which are defined in Ag^^, which is eo-close to the cylinder Aq 
obtained in Corollarv 13.21 

4.1 The specific form of the inner and outer maps 

For eo small enough the system associated to the Hamiltonian ()17p has a normally hyperbolic 
invariant cylinder Agg , which is eo-close to Ao given in Corollary 13. li Analogously, the Poincare 
map Veo associated to this system and the section {g = 0} has a normally hyperbolic invariant 
cylinder Ae^ = Aeo ^ {5 = 0}. Moreover, it is formed by seven connected components Agg, 
j = 0, . . . , 6, which are CQ-close to the cylinders Aig obtained in Corollarv 13.21 

Recall that in Corollary 13.21 we have seen that in the invariant planes / = constant there 
where forward and backward transveral heteroclinic connections between Aq and Aq provided 
I S and / G respectively. For the elliptic problem and eo small enough we will have 
transversal heteroclinic connections in slightly smaller domains. To this end, we define 

V*s = {I eV* : dist(I, dV*) >6}, * = f, b. (43) 

Theorem 4. For any 6 > 0, there exists eg > such that for < cq < e^ the map Vj^ , where 
Veo is the Poincare map associated to Hamiltonian (jlZp and section {g = 0}, has seven normally 
hyperbolic locall'!^ invariant manifolds A^o, which are e^-close in the sense to A^. Moreover, 
there exist functions Qi^ : [/_ + 6, 1+ — 6] x T — )• (M x T)"^, j = 0, . . . , 6 which can be expressed 
in coordinates as 

Qi,{i,t) = [gi;^ii,t),gi^ii,t),gifii,t),o,i,t) , (44) 

that parameterize Aig . In other words Ai^ is a graph over (/, t) defined as 

Ai^ = {^eo (/, t) : (/, t) G [/„ + 6,1+ -6] XT}. 

The invariant manifolds VF"(Ag^) and VF*(Ag^) intersect transversally provided I € 2?^ and 
the invariant manifolds W^{A'^^) and W''^{A^^) intersect transversally provided I G V^. More- 
over, one of these intersections is e^-close in the sense to the manifolds V^^^ defined in 
Corollary 13.^1 Denote T^eo these intersections. There exist functions 

c; (/, t) = (q^(/, t),c:^\i, t),c:f{i, t),o,i,ty * = f , b 

''see remark right below 
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which parameterize them. Namely, 

Ko = Koil^t) : {I,t) G [/_ + 5,/+ - 5] X T} , * = f,b. 

For the eUiptic problem, the coordinates (/, t) are symplectic not with respect to the canonical 
symplectic form dl A dt. Indeed, if we make the pullback of the canonical form dL Ad£ + dG A 
dg + dl f\ dt into the cylinders AgQ, we obtain the symplectic form 

= (l + eoai(-^, t) + ela^^{I, t) + e|]ai(/, t)) dl A dt, (45) 

for certain functions O'^ : [/_,/+] x T — )■ M. The functions a-L, contain the Bq remainder terms, 
and thus depend on cq even if we do not write explicitly this dependence to simplify notation. 
Remark 4.1. Objects and maps of Theorem H] have increasing regularity when eo tends to 0. 
Indeed, by Gronwall's inequality the Lyapunov exponents of Ae^ tend to zero with eo- So for 
every A; > 1, if eo is small enough, the invariant manifold Ae^ and the above subsequent objects 
are of class C'^ (see Remark 13. 6p . For the sake of simplicity, we will not henceforth emphasize 
regularity issues. But the main point is that for eo small enough all objects of our construction 
are smooth enough, and in particular it will be possible below to apply the KAM theorem to 
the invariant manifolds Aq. 

Remark 4.2. Theorem H] only guarantees local invariance for Agg. Namely, the boundary might 
not be invariant. Nevertheless, later in Section [5] we will show the existence of invariant tori in 
A|q, which will play the role of boundaries of A|(,. Thanks to these tori, one can choose Ae^ to 
be invariant. For that reason, from now on we will refer AgQ as a normally hyperbolic invariant 
manifold. 

We introduce the following notation. 

Notation 4.3. For any function f with 2TT-periodic dependence on t, M{f ) is the subset of the 
integers corresponding to which t-harmonics of f (which may depend on other variables) are 
non-zero. 

In the invariant cylinder A^^^ given in Theorem 2] one can define inner and outer maps as we 
have done in Aq for the circular problem. Next sections are devoted to the perturbative analysis 
of these maps. We state here their main result. 

Theorem 5. The normally hyperbolic invariant manifold given in Theorem^ of the Poincare 
map Vl^ in ([I2]) has associated inner and outer maps. Moreover 

• The inner map is of the form 

yi-. f I\^f I + eoAi{I,t) + elA2{I,t) + 0{el) \ 

■^eo-^t )^yt + Ml) + eMl.t) + elT2{I,t) + 0{el) J' ^^^^ 

and the functions A\, A2, 7i, and T2 satisfy 

M (Ai) = {±1}, M {A2) = {0, ±1, ±2} (47) 
M (Ti) = {±1}, M {T2) = {0, ±1, ±2}. (48) 

• The outer maps are of the form 
-po.t,*.( I + eoB*{I,t) + 0{el)\ 



eo 

and the functions B* satisfy 



M{B*) = {±1}. (50) 
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4.2 The Co-expansion of the elhptic Hamiltonian 

We obtain expansions of AHdi in ()17p with respect to cq. These expansions will be used later 
in Sections 14.31 14.41 and 14.51 The most important goal is to see which harmonics in t have the cq 
and Cg terms of the elliptic perturbation in ()17p . Note that the circular problem is independent 
oft. 

We define the function 

B(r,v, g,t) = -. — — — 7^^^ — — — — -— j-. (51) 

re*(^+^'-*) - ro(t)e™oW 

This function is the potential 

1 



k - Qo{t)\ 

expressed in terms of g = g — t, where g is the argument of the perihelion, the true anomaly 
V of the asteroid defined in ([13]) and the radius r. The functions ro(t) and vo(t) are the radius 
and the true anomaly of Jupiter. The functions ro(t) and vo{t) are the only ones involved in the 
definition of B which depend on cq ■ 

Then, the perturbation in (fTUj) can be expressed as 

1 — ^ j T 

nAHcnciL,i,G,g) + fieoAHcn{L,i,G, g,t) = B i — ,v,g,t 



^ -B (-^,v,g-t,t] + ^ 



{r,v)={r{L,e,G),v{L/,G)) 



1 — fi \l — fi J r 

First we study the expansion of B, and from it, we deduce the one of AiJeii- 
Lemma 4.4. Consider the expansion 

B{r, V, g, t) = Bo{r, v, g) + eoBi{r, v, g, t) + eg^2(r, g,t) + (eg) (52) 
of the function B defined in (j5ip . Then, 

• Bo satisfies N{Bq) = {0}. 

• Bi satisfies N{Bi) = {±1} and is given by 

Bi{r, V, g, t) = . r (2 cos t - 3r cos{v + g + t) + r cos{v + g - t)) , (53) 

2A-i{r,v,g) 

where 

1 /2 

A(r, V, (7) = (r^ + 1 — 2r cos{v + g)) 

• B2 satisfies M{B2) = {0,±1, ±2}. 

Note that the elliptic problem is a peculiar perturbation of the circular problem in the sense 
that the /c-th eo-order has t-harmonics up to order k. This fact will be crucial when we will 
compare the inner and outer dynamics in Section O 
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Proof. To prove this lemma we look for the eo-expansions of the functions ro(t) and vo{t) involved 
in the definition of B in (|51|) . We obtain them using the eccentric, true and mean anomalies of 
Jupiter. 

From the relation t = uq — cq sinuo (see (|15|) ). one can obtain that 



p2 



UQ{t) = t + eosint + sin2t + O (cq) . 
Then, using tq = 1 — eo cos uq, 

ro(t) = 1 - eo cos i + eg sin^ t + O (cq) . 
For the eccentric anomaly we can use 



vq / 1 + eo ^ uo 

tan — = \ tan — 

2 V 1 - eo 2 

(see ipS]) ). to obtain 

■^0 = ^0 + Co sin Uo + Cg sin uq — 2 sin 2iio^ + O (cg) 

and then 

VQ{t) = t + 2eo sint + eg Q sint - sin2t^ + O (eg) . 

Plugging rg(t) and V(){t) into ()5ip . it can be easily seen that the expansion ()52p satisfies all the 
properties of Bq, Bi and B2 stated in the lemma. □ 

From the expansion of B, one can easily deduce the expansion of in ([Ij 

Corollary 4.5. Let us consider the eo-expansion of the function AHpn in ()17p . 

AiJeii = AHl,, + eoAi/^i + O {el) . 

Then, the first order is given by 



AHl^,{L, i, G, g,t)=- ^Bi ( -^^hll^^y^L, i, G),g, t) 



Bi — ,v{L,i,G),g,t 



where Bi is the function defined in Lemma \4-4\ o-nd therefore it satisfies 

N[AHl,,)={±l}. 



(54) 



The second order A//^,, satisfies 



ell 

AA(A/72i)={0,±l,±2}. 
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4.3 Perturbative analysis of the flow 

To study the inner and outer maps perturbatively, we need first to study the first orders with 
respect to eo of the flow <&eo{s, {g, I, t)} associated to the vector field ([^T]) . Particularly, we want 
to know their dependence on the variable t. Recall that we already know the dependence on t 
of the 0-order thanks to formulas (1321) and (|3^ 



Lemma 4.6. The flow ^eo{s, {L,£,G, g, I,t)} has a perturbative expansion 

$eo{5, {L, £, G, g, /, t)} =$o{s, (L, G, g, I, t)} + eo$i{s, (L, £, G, g, I, t)} 

+ el'^2{s,iLJ,G,g,I,t)}+0{el), 

which satisfies 

Ar{^,{s,{L,i,G,g,I,t)}) ={±1} (55) 
AT {^2{s, {L, £, G, g, I, t)}) = {0, ±1, ±2}. (56) 

Proof. Let us denote z = (L, £, G, g, I) and call ^Ygy to the vector field (|41|) . which has expansion 

Xeo =Xo + eoXi + elX2 + O (ej]) . 

We first prove (j55p . The eo-order <I>i is a solution of the ordinary differential equation 

= DXo {<l>o{s, {z, t)}) e + Xi (Ms, {z, t)}) 

with initial condition ,^(0) = (0, 0). By ()3ip . Xq is independent of t and therefore, 

DXo {Ms, (z, t)}) = DXo ($r (s, z)) , 

where has been defined in ()32p . Then, this term is also independent of t. Prom Corollary 
I4.5[ one can deduce that M^Xi) = {±1} and therefore Xi can be written as 

Xi{z,t) = X+{z)e'^ + X-{z)e-'\ 

Therefore, using formulas (f32|) and (f37|) . one has that 

X, ($o{s, {z, t)}) = (X^ m'^'is, z}) e^*o{^'-}) e^* + (x^ (l>r{s, z}) e^*o{^'-}) e"^*. 



To prove ([55]) . it is enough to use variation of constants formula. Let us consider Mz{s) the 
fundamental matrix of the linear equation 

^J = DXo{<^>ris,z))^, 



then 
with 



(z, t)} = a>+{s, z}e'' + z}e-'' 
z} = MM M-\a) {Xt ($r{s, z}) e±^5oi^'^») da. 
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The proof of ()56p follows the same lines. Indeed, $2 is solution of an equation of the form 

= DXo ($r{5, z}) e + g, I, t) 

with initial condition ^(0) = (0, 0, 0). The function H is given in terms of the previous orders of 
Afeo and <I>eo as 

H = ]^D^Xq {^f') + DXi {^f') + X2 > 

and then it satisfies M{E) = {0, ±1, ±2}. 

Since the homogeneous linear equation is the same as for $™ and does not depend on t, one 
can easily see that (j56|) is satisfied. □ 



4.4 Perturbative analysis of an invariant cylinder and its inner map 

We devote this section to study the normally hyperbolic invariant manifold of the elliptic problem 
Agjj, given in Theorem U and the inner map associated to it. We study the inner map of the 
elliptic problem perturbatively from ()26|) . taking cq as a small parameter. We call the inner 
map J-"g" : A^^^ — )• Ag^, which is defined as the (— 147r)-Poincare map of the flow ^>gg, deflned in 
Lemma restricted to Ag^^. 

We want to see which i-harmonics the first orders of the inner map has and also we want 
to compute the first order of the /-component. To this end we consider the classical theory of 
normally hyperbolic invariant manifolds |Fen741 IFen77j . This theory ensures the existence of 
the maps parameterizing the normally hyperbolic manifolds Aq of the map 'Pj^^ . Moreover, 
they can be made unique imposing 

7rigi^{I,t) = I,7rtgi^iI,t) = t, (57) 

where n^, is the projection with respect to the corresponding component of the function. Since 
we only need properties of the parameterization of Ag^ and the dynamics on it, we just consider 
the case j = 3. The map Q^^^ satisfies the invariance equation 

reoogi = gloTr^, (58) 

where Veo = "Pj^ and J-"™ is the inner map of the elliptic problem, namely the Poincare map Vj^ 
restricted to the cylinder Ag^^. 

Since we have regularity with respect to parameters, the invariance equation allows us to 
obtain expansions of both the parameterizations of Ag^ and the inner map J-™ with respect to 
Co- Let us expand and .F™ as 

K = ^^ + eo-^r + elTt + O {el) . (60) 

Then, is the function defined in (j24p and J-"™ is the inner map of the circular problem obtained 
in ()26p . which is defined in Ajj. Recall that 

Ve,{LJ,G,0,I,t)=rLiL,£,G,0,I,t)='^eo{-l^^,{L,i,G,0,I,t)}, (61) 
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where is the flow considered in Lemma 14.61 Then, we have that 

M (Vi) = {±1} and N {V2) = {0, ±1, ±2}. 



Then, expanding equation (j58|) with respect to eo allows us to deduce the properties we need 
about the inner map. They are summarized in the next lemma, which reproduces the part of 
Theorem [5] referred to the inner dynamics. 

Lemma 4.7. The expansions of the functions Q'^^ and J^g" in ()59p and ()60p satisfy that 

N{gl)={±l} and (gf ) = {0, ±1, ±2} 

and 

AA(^{") = {±1} and M{F'^) = {0, ±1, ±2}. 
Namely, the inner map is of the form, 

^\^( ^ + ^0^1 ^) + e^^2 (/, t)+0 {^el) \ 
■^eo-^t y ^ ^ i + ^7-Q(j) + eQri(/,i)+e2r2(/,t) + 0(/ie2) )^ ^'''> 

and the functions A\, Ai T\ and Ti satisfy 

M{Ai) = {±1}, {A2) = {0, ±1, ±2} (63) 
M{Ti) = {±1}, M {T2) = {0, ±1, ±2}. (64) 

Moreover Ai can he split as, 

A,{I,t)=AtiI)e'' + A^iI)e-^\ 

with 

AtH) = .t^'FAH^ti . ^^^Ua, (65) 

Jq -1 + ^(9GAi?circA|(cr) 

where the functions /S.H^^^ are defined as 

AHl,,iL, i, G, g, t) = AHl^iL, i, G, g)e'' + AH^iL, £, G, g)e-'\ 
A^((t) has been defined in ([33j) and 

\]{a) = l,{a,{Gl'\l),gl\l),gl^''{I),Q)}, (66) 

where has been defined in ()72p and Qq in Corollary \3.S\ . 

From the properties of ^f^, we can deduce properties of the symplectic form Og^^ in (fl5]) . 
which is defined in the cylinder Ag^^. 

Corollary 4.8. The functions a\ and in the expansion of the symplectic form in ()45p . 
which is the puUback of the symplectic form dL A d£ + dG A dg + dl A dt into the cylinder A^^ , 
satisfy 

M{al)={±l} and M {afj = {0,±l, ±2}. 
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Proof of Lemma \4- ?[ In the proof we omit the superscript 3 of the terms in the expansion of 
Expanding equation ()58|) with respect to cq, we have that the first terms satisfy 



Vo o Go =go o (67) 
Pi o ^0 + (dVo o ^o) Gi =Qi o K + [DQo o K) (68) 



+ ( DPO o go ) ^2 =^2 o + [DGi o Jj-) ^> 



2 

By the uniqueness condition (|57|) . t/i is of the form 



+ i (Z^^^o o J-^") (^r)®2+ (i^^o o J-i-. (69) 



GiigJ^t) = (^gi(5,/,t), 0,0,0 

with gi{g,I,t) = {g^{g,I,t),gi{gJ,t),gf{g,I,t)). 

Equation ()67p corresponds to the inner dynamics of the circular problem. We use equations 
(|68p and (|69|) to deduce the properties of J^™ and J^™ respectively. These equations can be 
solved iteratively and thus we start with ([68|) . Since 

ngo = ) and Z)a. = ( ) for i > 1, (70) 

we have that 



= ^* [Pi ° ^0 + [DVo o ^0 j ^1 j , * = /, t. 

Replacing this into (|68p we obtain an equation for t/i . The equation for every Fourier ^-coefficient 
is uncoupled, and therefore, using ([6T]l . that Vq is independent of t and the uniqueness of ^i, 
one can deduce that Mi^gi) = {±1}- As a consequence we have that AA(J^"^) = {±1}. 

Reasoning analogously and using ([6T]) again, one can see that M{g2) = {0, ±1, ±2} and 
Ar(J-^°) = {0,±l,±2}. 

Now it only remains to prove formula (|65p . To this end, let us recall that one can write the 
/-component of the inner equation as 

J-lf (/,t)=<{-147r,geoU,t)} 

since it is defined as the (— 147r)-Poincare map associated to the fiow of system (j4ip restricted to 
the cylinder Ag^ . Recall that the minus sign in the time appears due to the fact that system ()4ip 
has the time reversed with respect to the original one. Then, one can apply the Fundamental 
Theorem of Calculus and use ()4ip to obtain 

^147r 



(/,i) = -<{.,geo(/,t)}ds 



/ieo9fAgell0^eo{^,geo(/,t)} 

-1 + ^idcAHcirc O $eo {S, geoil, t)} + lieodgAH^n O ^eo {s, geo {I, t)} 



ds 
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Then, using the expansions of the Hamiltonian Ai^eii in CoroUarv l4.5t of the flow $eo ™ Lemma 
14.61 and of Qeo just obtained, 

- Jo -l + /.acAF_o$o{.,ao(/,i)}^' + ^ 

Namely, 

' ^ io -l + A^9GA//circO<i>o{s,ao(/,i)} 
To deduce the formulas for Af it is enough to split Ai^^jj as 

AH^,,{L,i,G,g,t) = AHl^iL,i,G,g)e'' + A/7^\f (L, ^, G, g)e-**, 
and recall that by formulas (|32p and ()37p . <l>o can be written as 

$0 {s, {L, e, G, g, I, t)} = (<i>circ {s, {L, £, G, g, I)},t + $o {s, (L, ^ G, g, I)} ^ 

□ 



4.5 The outer map 

We devote this section to study the outer maps 

-^^'^^^eo^^eo, * = f,b (71) 

for eo > 0. 

Theorem m gives the existence of C*^, * = f,b, transversal intersections of the invariant 
manifolds of Ag^^ and A^^. The Poincare map ()42p possesses then also a normally hyperbolic 
invariant manifold Agp = U^^qA^q and its invariant manifolds intersect transversally at Te^. 

Then, we can proceed as in Section 13.31 to define the outer maps T^^K We want to study 
it as a perturbation of the outer map of the circular problem given in ()30p . To this end we use 
Poincare-Melnikov techniques. As we have explained in Section [3.31 the original flow associated 
to Hamiltonian (fTBj) does not allow us to study perturbatively J^°o^- Instead, we use the reduced 
elliptic problem defined in (|4ip . 

The results stated in Theorem [5] about the outer map are summarized in the next lemma. We 
also show how to compute the first order term. To this end, we use the notation used in Section 
13.31 In particular, we will use the notation 7^'^(o") and Xf^{a) defined in ([33]) . Analogously we 
define their corresponding t-component of the flow as 

7K^) = 5o{^, (Co*'^/),C*'^(/),Co*'^(/),0)}, * = f,b 

x^ia) = $o{^, ig^'\i),g^'\i),gi'''ii),o)}, j = 3,4 

where has been defined in (j37p and Cg and Gq have been given in Corollarv l3.2[ 
Lemma 4.9. The outer map defined in ()7ip has the following expansion with respect to cq, 



(73) 
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Moreover, the functions B*'^{I) can be defined as 

B'^'^il) = B^{I) + B^^i{I)e^'^<('\ 
where uj^^^^^{I) and ujfj^il) are the functions defined in (|34p and (|36p respectively and 

Bl;^,{I)=±^^^ lim N ^^Ij o ^\{a) 

-l + ^^cAi^circO A3((j) J 



T^/. lim / ^ifoiLJiin^^ginlH 



-l + //c)GAi/eircOAf(cT) i ' 



where 



AHl^ii, L, g, G, t) = AHl^ii, L, g, G)e'' + A/7^\f (£, L, 5, G)e- 
/las &een defined in Corollary |^.5| antf Wj- /laue 6een defined in (j35p . 



(74) 



(75) 



5„^^(/) = ± z/. lim r f ^^^K^^^^^Kr-.^^'''^''^ 

^^elf °^K^) g±i(A|M+M<(/))\ ^ .^g^ 

-l + /.aGAi/eircoAKa)' J'''' ^'^^ 



(77) 



Proof. To prove the lemma recall that the outer maps are the composition of two maps. Indeed, 
as we have explained in Section 13.31 they are defined as 

•r-out.f _ ^6 (-f .73 , T3 

eo — ' eo ° '^eo ' ^^0 ^ ^^0 

•770ut,b _ ch . T3 , 73 



Thus, we will study both maps perturbatively and then the composition of them will lead to the 

r-OUt,f ■ ; 1 r r T-OW 

'"eo smce the prooi for y-go 



proof of the lemma. We only deal with J^"*'^ since the proof for J-^^*'^ can be done analogously. 
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To study Sl^ : Aq — Aq we use the Definition 13.41 of the (heteroclinic) outer map. Let us 



consider points z G r*^, £ A^^^ and x_ G A^^^ such that 



dist {V2,(,z),V2,{x±)) < C7A^I"I for n G Z± 

for certain constants C > and A > 1. Using the parameterizations of Fg^ and Ai^, j = 3,4, 
given in Theorem HI we can write the points z and x± in coordinates as z = Ceo(/o,to)) = 
Q^^{I+,t+) and X- = G^^{I-,t-). Then, the /-component of the outer map is just given by the 
relation 

= /+ = /_ + (/+-/_). 

To measure /+ — /_ we first deal with Iq — I±. Consider the flow associated to the reduced 
elliptic problem (|4ip . Then, applying Fundamental Theorem of Calculus 



/0-/4 
lo-I- 



Note that the changed signs in the limits come from the fact that system ()4ip has the time 
reversed. 

Using the perturbative expansions of C^^ and Agg given in Theorem [H equation ()4ip . the 
perturbative expansion of the Hamiltonian (jl6p given in Corollary 14.51 and the perturbation of 
the flow given in Lemma l4.6t one can see that 



/n 



-eo lim 



lidt^Hl,,{L,l,G,g,t) 



T 



-l + fidGAH,i,,{L,i,G,g) 
l^dtAH^,,{L,£,G,g,t) 



(L,£,G,g,t)=(<S>-^'^-,'fl){s,CliIo,to)} 



-eo lim 

T-!>+oo _/y 



-l + fidGAH,i,,{L,i,G,g) 
V fxdtAHl,,{L,e,G,g,t) 



-l + fidGAH,i,,{L,i,G,g) 
fidtAH^,,{L,£,G,g,t) 



{L,e,G,g,t)={%^^'=,<S>i^){s,gl{I+,t+)} J 
(L,£,G,g,t)=('S>-^'^-,'fl){s,CliIo,to)} 



ds + O [el] 



-l + fidGAH,i,,{L,i,G,g) 



{L/,G,g,i)=(<I.g'-,**){.,giJ(/-,t-)} 



\ds + {el) 



where <I>q ''^ and $q have been defined in ()32p and ()37p respectively. 

Taking into account that in Corollarv 14.51 we have proved that AH^^^ satisfies A/" (Ai/^^jJ 



{±1}, one can easily obtain the formula for 5^1^ in ([75]) . 



To obtain the formula for one can proceed as in the study of the inner map in Section 
14.11 Finally, to obtain the formula for B^'^ it is enough to compose both maps and . □ 



5 Existence of diffusing orbits 

The main result of this section is the following. 

Theorem 6. For any 6 > there exists Bq > and C > such that for any < eo < Cq the 
map Veo in (|42p has a collection of invariant 1- dimensional tori {Ti}^^ C Ag^ such that 
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• Ti n {/ = /_ + 5} / and T^v n {/ = /+- 5} / 0. 

• Hausdorff dist{Ti,Ti^i) < Ccq^^. 

• These tori form a transition chain. Namely, Wj, ftl 7^ for each i = 1, . . . , iV — 1. 

Proof. Once we have computed the first orders in eo of botli the outer and inner maps, we want 
to understand their properties and compare their dynamics. To make this comparison we will 
perform two steps of averaging |AKN88[ ILM88j . This change of coordinates will straighten the I- 
component of the inner map at order ©(bq) in such a way that with the new system of coordinates 
will be easier to compare the dynamics of both maps. Nevertheless, before performing averaging, 
we have to perform a preliminary change of coordinates to straighten the symplectic form 
to deal with the canonical form dl A dt. 

Lemma 5.1. There exists a e^-close to the identity change of variables 

{I,t) = {l',t')+eo^i{l',t'), (78) 
defined on A^^, which transforms the symplectic form defined in (02]) into the canonical form 

Qo = dl' A dt'. 

In the new coordinates, 

• The inner map J-"™ in ()46p reads 



I' + eo^i (/', t') + elA^ (/', t!) + O [lief) 
t' )^\t'^ ^To (/') + e^T{ {I',t') + elT^ (/', t') + O (/xe^) 



where A\ is the function given in Lemma \J^.l\ and N^, T[ and 7^' satisfy 

AT {A'2) = {0,±1,±2}, Ar{Tl)={±l} and AT {T2) ={0,±1, ±2}. 

• The outer maps J-^o^'^ and Je*"*'^ in ()i9]) read 



where B* are the functions given in Lemma \4-9[ 
Proof. We will see that there exists a change of coordinates of the form 

(i = r + eif2{r,t') + o{ef,) 

[t = t' + eogi (/', t') + elg2 (/', t') + O (e^) 



(81) 



with 

M{gi) = {±1}, M{g2) = {0, ±1, ±2} and M{f2) = {0, ±1, ±2}, (82) 

which straightens the symplectic form rig^^. In fact, we look for the inverse change. Namely, we 
look for a change of variables of the form 

{r = I + elh{I,t) + elJ>{I,t) 
\ t' = t + eogi{ht) 
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such that the pullback of Oq = dl' A dt' with respect to this change is the symplectic form ilg^^ . 
Even if we do not write it expHcitly, /> depends on cq. To obtain this change, it is enough to 
solve the equations 

dtgi = al 
dih = 4 
dif> = b 

where 

6 = a| - dtgidif2 - eodtgidif> + digidtf2 + eodigidtf> 

and af , and a> are the functions considered in ()45p . These equations can be solved iteratively. 

Recall that by Corollarv 14.81 we have that A/'(af) = {±1}- Then, one can take gi as the 
primitive of af with zero average, which satisfies 

M{gi) = {±l}. (84) 

The other equations can be solved taking 



/2(/,t)= f al{J,t)dJ 
Jo 

7>(/,t) = / b{j,t)dj. 

Jo 



Note that b depends on gi and /2, which have been already obtained. Since by Corollarv 14.81 we 
have that AA(a|) = {0, ±1, ±2}, one can deduce that 

ATf/a) ={0,±1,±2}. (85) 



To obtain the change ()8ip it is enough to invert the change ()83p . Then, formulas ()84p and ()85p 
imply (HI]). 

To finish the proof it only remains to check the properties of the inner and outer maps in 
the new coordinates. To this end one only needs to take into account ()82p . □ 



Once we have straightened the symplectic form, we can perform two steps of averaging to 
the inner map. 

Lemma 5.2. There exists a CQ-close to the identity symplectic change of variables 

{l',t') = {I,r) + eoip2(I,r) (86) 

defined on A^^, which: 

• Transforms the inner map 

J^'eo into 

■ U y U + mTo (X) + elT2 (X) + O {^,el) J ^''^ 



36 



Transforms the outer maps ' and Fea ' in ([80]) into 

Tont,* .( ^\^( ^ + eoB* {I, t) + if^el) \ 
V r + ^a;*(X) + 0(/ieo) J' 



where 



with 



B* (I, t) = B*'+ (X) e'^ + S*'^ (X) e~'^ 

= _ -^-^^^^Af ix) . 

Moreover, these functions satisfy 

B*'^{I)^0 for lev*, (89) 
where T>* are the domains considered in Corollary[ 



Note that we can do two steps of averaging globally in the whole cylinder AgQ due to the 
absence of resonances in the first orders in eo- Namely, there are no big gaps. This is quite in 
contrast with the typical situation in Arnold diffusion (see e.g. jDdlLS06] ). 

Proof. We perform two steps of (symplectic) averaging. To this end we consider a generating 
function of the form 

S{I, t') = It' + eo5i(X, t') + eg52(X, t'), 
which defines the change ()86p implicitly as 

1 = 1 + eo9t'5i(X, t') + eldt,S2{I, t') 
T = t' + eodzSiil, t') + eldxS2{T, t'). 

By Theorem 13. 3l we have ()27p and by Theorem [5] we know which t'-harmonics have the functions 
Ai and %■ Then, it is easy to see that the functions Si that correspond to two steps of averaging 
are globally defined in A^^^. Moreover, in these new variables, taking into account that the inner 
map is exact symplectic, one can see that the inner map is of the form (j87p . 

To obtain a perturbative expression for the outer maps J-?o"*'*, we need to know Si explicitly. 
It is given by 

9 rx n = _J^l^,^t , ^-47 (X) 

If we apply this change to the outer maps 1^^^'* in ([73]). we obtain ([88|) . 

The statement (|89p is based on convincing numerical data (see Appendix IC.3p . □ 

In the new coordinates (X, r) the inner map in ()87p is a eg-close to integrable map. 
Moreover, thanks to Theorem 13.31 is twist and therefore we can apply KAM theory to proof the 
existence of invariant curves in Ag^^ . We use a version of the KAM Theorem from [DdlLSOO] (see 
also |Her83p . 

KAM theorem. Let f : [0, 1] x T — t- [0, 1] x T he an exact symplectic map with I > 4. Assume 
that f = fo + 6fi, where fo{I,i;) = + A{I)), A is C\ \diA\ >M and\\fi\\ci < 1. Then, if 
8^I'^M~^ = p is sufficiently small, for a set of oj of Diophantine numbers of exponent = 5/4, 
we can find invariant tori which are the graph of C^^^ functions u^, the motion on them is C'~^ 
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conjugate to the rotation by ui, and \\ui_j\\ci--i < C(5^/^. 

Applying this theorem to the map we obtain the KAM tori (see Remark 14.11 for the 
matter of regularity). Moreover, this theorem ensures that the distance between these tori is no 
bigger than eg^^. Then, the results of Lemma 15.21 and KAM theorem lead to the existence of a 
transition chain of invariant tori. 

The transition chain is obtained comparing the outer and inner dynamics. We do this 
comparison in the coordinates (X, r) given by Lemma 15.21 and therefore we deal with the maps 
Tg" and J-?o^*'* in ([87|) and (f88|) respectively. 

KAM theorem ensures that there exists a torus Ti such that Ti n {/ = /_ — 5} 7^ 0. Then, 
either J-^"*'^ or 7-"°"*''^ are defined for points in Ti. Assume, without loss of generality that J-^q^'^ 
is defined for points in Ti. Then, thanks to (|89p . J-"e"*'^(Ti) satisfies 

-out,f 



dist (^Ti,J-°;^*''(Ti)j >Ceo 

for a constant C > independent of cq. Then, KAM theorem ensures that there exists a torus 
T2 such that T2 r\ F^^^'^ {Tij / 0. Iterating this procedure, choosing at each step either Tt"^'^ 
or J-^"*'*^, one obtains the transition chain. □ 

To finish the proof of Theorem [2] it only remains to prove the existence of a diffusing orbit 
using a shadowing method. We use a result stated in [FMOOj . 

Lemma 5.3. Let f be a sympletic map in a 2{d + 1) syniplectic manifold. Assume that the 
map leaves invariant a d- dimensional torus r and the motion on the torus is an irrational 
rotation. Let T be a d + 1 manifold intersecting transver sally. Then, 



W^c[jf-^iT). 



i>0 

An immediate consequence of this lemma is that any finite transtion chain can be shadowed 
by a true orbit, which finishes the proof of Theorem [2j 

In terms of the elliptical elements of the asteroid, such a diffusing orbit can now be described 
as follows. The orbit starts near the resonant cylinder A. The eccentricity of the primaries is 
small: this is an essential feature of both the proof above and the qualitative behavior of the 
orbit. So, over a time interval of length <^ I/cq, the orbit closely follows a hyperbolic periodic 
orbit of the circular problem. The semi major axis is roughly constant equal to 7^^^ and the 
Jacobi constant to —1.81. The asteroid turns around the primaries, making one full turn over a 
time interval of length 7. In the frame rotating approximately with the primaries, the Keplerian 
ellipse itself of the asteroid rotates counterclockwise with fast frequency ~ —1; in the inertial 
frame of reference, it rotates only /x-slowly (see e.g. |AKN88l Fej02| ), while the eccentricity e 
slowly oscillates around e = 0.48. 

At some point (as soon as we can if we want to save time), the orbit undergoes a hete- 
roclinic excursion, during which a heteroclinic orbit is shadowed over a time interval of size 
0(— ln(/ieo)/y7^). During this excursion, the semi major axis itself undergoes an oscillation of 
magnitude 0{^/Jl), eventually coming back to its initial approximate value 7^/^. On the other 
hand, the Jacobi constant and the eccentricity have increased by O{^eo). Excursions may then 
repeat, piling up changes in the eccentricity to reach the value e = 0.66 in finite time. 
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A Numerical study of the hyperbolic cylinder of the circular 
problem. 

We devote this appendix to the numerical study of the hyperbohc invariant manifold of the 
circular problem given in Corollary 13.11 and its invariant manifolds. In other words, we show 
numerical results which justify the properties of the circular problem stated in Theorem [3l 

Numerical analysis has several sources of error: mainly round-off errors in computer arith- 
metics, and approximations of ideal mathematical objects (e.g. linear approximation of local 
stable/unstable manifolds). In our analysis, we have tried to evaluate such errors, and check 
that they are appropriately small. We do not claim to give a fully rigorous proof of Theorem [3l 
which would require Computer- Assisted techniques as in |WZ03j . Indeed, we have focused our 
efforts to keep the numerics relatively simple and, hopefully, convincing. One could think of 
several possible numerical computations to prove our result. The most numerically demanding 
one would be to check directly that some given orbit has an adequate drift of eccentricity. This 
computation would not bring much light to the mechanism of instability, and moreover it would 
involve formidable numerical analysis problems, due to the necessarily very long time of integra- 
tion. On the contrary, our line of proof allows us to use numerical verifications involving only 
orbits of the circular problem -a dramatic simplification, as we will see below. 

Let us make a few more specific comments on the strategy of our numerical analysis. As 
mentioned in Section O the circular problem has a conserved quantity, the Jacobi constant which 
we denote by J (see ([5])), which corresponds to energy when the system is expressed in rotating 
coordinates. Thus it is natural to fix the Jacobi constant J = Jq and perform our analysis for 
a given Jq. This allows us to reduce the dimension of the computations by one. Finally, we let 
J vary and repeat the computations for all J in the range of interest J S [J_, J+]. 

Another important comment is on the choice of coordinates. For numerics we prefer Cartesian 
coordinates, since the equations of motion are explicit in these coordinates. Thus we carry out 
our computations of the hyperbolic structure of the circular problem in Cartesian ( Appendix |A|) . 

On the other hand, for perturbative analysis we have used Delaunay coordinates throughout 
this paper. Thus, in Appendix [B] we explain how to change coordinates from Cartesian to 
Delaunay, and we carry out our computations of the inner and outer maps of the circular and 
elliptic problems in Delaunay (Appendix lUj). 

Regarding the integration method, we use a variable-order Taylor method specially gener- 
ated to integrate the equations of motion and variational equations of the circular problem. 
The Taylor method has been generated using the "taylor" package of A. Jorba and M. Zou 
(see http: //www.maia.ub. es/~angel/taylor/p . The main advantage of using a Taylor method 



is that it is very fast for long-time integrations (without sacrificing accuracy). 

A.l Computation of the periodic orbits 

Consider the circular problem in rotating Cartesian coordinates 

J{x,y,Px,Py) = ^(pI +pl) + ypx-xpy - - A (90) 



where 



2 f \2 I 2 

2 r , \2 I 2 

r2 = {x + fii) +y . 
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Recall that the energy of the circular problem in rotating coordinates coincides with the Jacobi 
constant J in ([S]). From now on in this appendix we will refer to J as the energy of the system. 

We follow the convention to place the large mass (Sun) to the left of the origin, and the 
small mass (Jupiter) to the right. (This is opposite to the astrodynamics convention). Thus we 
choose /ii = /i as the small mass, and ij,2 = 1 — n as the large mass. Notice that equation ([90]) 
is reversible with respect to the involution 

R{x,Px,y,Py) = {x,-Px,-y,Py)- (91) 

Thus, a solution of the system is symmetric if and only if it intersects the symmetry axis 
Fix(i?) = {y = 0, Px = 0}. This symmetry will facilitate our numerical computations. Note 
that the involution R is just the involution (j20p expressed in rotating Cartesian coordinates. 

Let the energy be fixed to J = Jq. We look for a resonant periodic orbit Ajq of (|90p 
in the level of energy Jq. As a first approximation to Aj^, we look for a resonant periodic 
orbit of the two-body problem, i.e. of the Hamiltonian (|19p with /_i = 0. Let us denote the 
approximate periodic orbit by A jg = {L,£,G, g). The actions L and G are determined by the 
resonant condition = 7 and energy condition — ~ G = Jq. To determine Aj,j completely, 
we choose that the asteroid is initially at the perihelion, i.e. we impose an initial condition 
Ajg(O) = {L^ ,1^ .,G^ ,g^) with = and = 0. Switching to Cartesian coordinates, we obtain 
an initial condition (x'^,p|^, y'^,Py) with = and y*^ = 0. 

Next we refine the trajectory Ajq into a true periodic orbit Ajq for the system ()90p with 
^ = 10"'^. Consider the Poincare section 

S+ = {y = 0, py > 0} 

in the circular problem (|90p . and let P : S"*" — )• S"*" be the associated Poincare map. Since we are 
in rotating coordinates, this section corresponds to collinear configurations of the three bodies. 

Remark A.l. In numerical integrations, we use a variable-order Taylor method with local error 
tolerance 10~^^. Moreover, a point is considered to be on the Poincare section whenever \y\ < 
W-^^ and Py > 0. 

Furthermore, the momentum variable py can be eliminated. Indeed, since dpyJ 7^ 0, py in 
the region of the phase space we deal with, it can be recovered from the other variables using the 
energy condition J{x,px,y,Py) = Jq. Hence, the Poincare map is a two-dimensional symplectic 
map at each energy level, acting only on {x,px). 

Notice that, in the rotating frame, a 7:1 resonant periodic orbit makes 6 turns around 
the origin. See Figure [H In principle, we could look for the periodic orbit as a periodic point 
p = {x,px) of the Poincare map: p = P^{p). This would imply solving a system of two equations. 
Thanks to the reversibility (j9ip . in fact it is only necessary to solve one equation. Notice that 
our initial condition {x,px) is at the symmetry section {y = 0, Px = 0}, so the periodic orbit 
must be symmetric. Thus it is enough to impose the condition that the trajectory Ajo(i) after 
half the period is again at the symmetry section. Hence we set up the problem as simple one- 
dimensional root finding: we look for a point p = (x,0), such that its third iterate P^{p) has 
momentum px = 0: 

7rpAPHp))=0. 

(Here, vr^^ : — )■ M is the projection onto the Px component). 

In order to solve this problem, we use a a Newton-like method. Specifically, we use a modified 



version of Powell's Hybrid method (see the GSL manual G"*" for details) without scaling. In our 
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X 

Figure 5: Resonant periodic orbit A_i.6 of the circular problem in rotating Cartesian coordinates. 

computations, the Newton method converges in less than 5 iterations. As a test of the software, 
we have checked that the rate of convergence of the Newton method is quadratic. 

Remark A. 2. We ask for an accuracy of 10~^^ in the Newton method, i.e. a point p = (x,0) is 
accepted as a root if and only if its third iterate P^{p) has momentum \px\ < 10~^^. 

For the Newton method, we need to compute the derivative of the Poincare map. For each 
^ S M^, let u{t, be the solution of the system with initial condition u{0, ^) = ^. Let T : S"*" — )■ M 
be the Poincare return time. The derivative of the Poincare map at a point p S M"^ is given by 
the partial derivative DP{p) = u^{T{p),p). It is well-known that u^{t,p) is the matrix solution 
of the variational equation 

W = Df{u{t,p))W, 

where / is the vector field of the circular problem. We compute DP[p) by numerically integrating 
the variational equation using the Taylor method mentioned above. 

For illustration, let us show some numerical results corresponding to the energy value J = 
— 1.6. The first approximation A_i.6 from the two-body problem has initial condition p^ = 
{x^ ,p^) = (1.30253- •• ,0). After refining this initial condition via the Newton method, we 
obtain a resonant periodic orbit A_i.6 of the circular problem passing through the point p = 
{x,p^) = (1.29858- •• ,0), with period T.i.g = 44.01796--- ~ 147r. See Figure E The periodic 
orbit A_i.6 is symmetric, with the points p and P^{p) located at the symmetry section (they 
have y = and px = 0). Notice that, in rotating coordinates, the trajectory of the asteroid 
makes 6 turns around the origin before closing up at the point p. 

Finally, we let J change and, using this procedure, we are able to obtain the resonant periodic 
orbit for energy levels 

J G [J-,J+] = [-2.04, -1.56]. (92) 
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Figure 6: Resonant family of periodic orbits. We show normalized period Th — 147r, and 
maximum deviation of L component with respect to the resonant value 7^^^ (see equation ()93p ). 



See Figure [6l This family of resonant periodic orbits constitutes the normally hyperbolic invari- 
ant manifold Aq given in Corollary 13.11 Notice that the period Tj stays close to the resonant 
period 147r of the unperturbed system. From Figure El we obtain the bound 

\Tj - U7t\ < 60fi, 

which is the first bound given in Theorem [3l 

To determine the stability of the periodic orbit Aj^,, we compute the eigenvalues A and A~^ 
of the matrix DP^{p), where DP^{p) is the linearization of the iterated Poincare map about 
the fixed point p. 

Figure [7] shows the characteristic exponents ln(A), ln(A~^) as a function of energy. The 
family of periodic orbits is strongly hyperbolic as J — )• J+, and weakly hyperbolic as J — )• J_. 
Note that one would expect that we are in a nearly integrable regime since is small. Then one 
would expect the eigenvalues to be close to 1. Nevertheless, in this problem the non-integrability 
is very noticeable when one increases fj, to fi = 10~^. This is due to the effect of the perturbing 
body (Jupiter) on the asteroid, as the asteroid passes close to it. 

Furthermore, we verify that (the square of) the semi-major axis L stays close to the resonant 
value 7"^/^. Integrating the periodic orbit in Delaunay coordinates Xj{t) = {Lj{t),ij{t), Gj{t),gj{t)) 
over one period Tj, we compute the quantity 

imax(J)= max \Lj{t)-7^/^\. (93) 

te[o,Tj) 

The function L^g^x{J) is plotted in Figure O Notice that we obtain the bound 

\Lj{t) - 7V3| < 7^ 
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Figure 7: Characteristic exponent ln(A) as a function of energy level J (the other exponent is 
-ln(A)). 
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Figure 8: Extremal periodic orbits of the family: circular periodic orbit with J = (in red), 
elliptical periodic orbit with J = (in green). The Lagrange equilibrium point L2 is marked 
with a '+' symbol. 
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for all t G M and J € [J_, J+J, which is the second bound given in Theorem [3l 

Let us briefly describe the family of periodic orbits Aj. For illustration, see Figure [51 At one 
endpoint of the family, as J — )■ J_, the periodic orbit Aj tends to a circular orbit of period 147r 
centered at the origin and passing far away from the primaries (Sun and Jupiter). Moreover, 
A J looses hyperbolicity when J — J„. For instance, the periodic orbit Aj_ of the two-body 
problem approximation has eccentricity e( J_) = 0.09989 • • • . 

At the other endpoint of the family, as J — )■ J+, the periodic orbit Aj tends to a homoclinic 
loop of the Lagrangian equilibrium point L2 that makes 6 turns around the Sun- Jupiter system. 
(In rotating Cartesian coordinates, L2 is located on the x axis at the point X2 — 1.068). This 
explains the fact that the period Tj "explodes" as J — >■ J+. Since we are interested in working 
close to the resonance, we avoid energies J > J+ where the period explodes. 

A. 2 Computation of invariant manifolds 

In this appendix, we compute the stable and unstable invariant manifolds associated to the 
periodic orbits found in the previous section. 

Consider first a fixed energy level J = Jq. Let Ajq be the resonant periodic orbit of the 
circular problem found in the previous section. To compute the invariant manifolds of the 
periodic orbit, we continue using the iterated Poincare map. Thus we look for (one dimensional) 
invariant manifolds of a hyperbolic fixed point at each energy level. Let p S Ajq be a hyperbolic 
fixed point of the iterated Poincare map P = P^. Let A, A"^ be the eigenvalues of DP{z) with 
A > 1, and Vu,Vs be the associated eigenvectors. 

Assume that we want to compute the unstable manifold IF"(p). Let 77 be a small displace- 
ment in the unstable direction Vu- We approximate a piece of the local manifold by the linear 
segment between the points p+rjVu and P{p+rjVu)- We call this segment a fundamental domain. 
We discretize the fundamental domain into an array of points, and iterate them by P to globalize 
the manifold. (The stable manifold is computed analogously using the inverse map P~^.) 

The error commited in the local approximation P{p+rjVu) = p+\riVu-\-0{rj^) of the manifold 
is given by 



Remark A. 3. For each energy level J, we choose a displacement r] = rj{J) such that the local 
error is err(?7) < 10~^^. 

One can think of p as a fixed point of the iterated Poincare map P = P^, or as a 6-periodic 
point of the Poincare map P. If pi = P^{p) are the iterates of p for z = 0, . . . , 5, then pi are also 
fixed points of P. They have associated unstable and stable manifolds, which can be obtained 
from W^'^{p) by iteration. 

For illustration, let us show some numerical results corresponding to the energy value J = 
— 1.6. Figure [9] shows the manifolds of all iterates {pj}i=o,...,5- Notice that the dynamics in 
Figure is reversible with respect to the symmetry section {y = 0, px = 0}, as discussed in 
the previous section (see ()9ip ) . Figure [9] shows that the manifolds do intersect transversally 
at different homoclinic points. We are interested in measuring the splitting angle between the 
manifolds. Unfortunately, the homoclinic points do not lie on the symmetry axis, which would 
be very useful in order to compute them. 

In order to have the homoclinic points lie on the symmetry axis, we recompute the manifolds 
on the new Poincare section 




S- = {y = 0, py < 0}. 
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Figure 9: Invariant manifolds of the fixed points PoiPii ■ ■ ■ :P5 on the section S"*". Unstable 
manifolds are plotted in red, stable in blue. The fixed points are marked in green. 
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Figure 10: Invariant manifolds on the section S~. 



45 




Figure 11: Invariant manifolds of the points p2 and on the section S . Due to the symmetry, 
points that he on the hue = (marked in green) are intersection points. 



Numericahy, we just transport points on the unstable manifold from section S"*" to section S~ 
by the forward flow, and points in the stable manifold by the backward flow. See Figures [10] 
and [TTJ Now the points that lie on the symmetry line = are homoclinic points. 

A. 3 Computation of transversal homoclinic points and splitting angle 

In this appendix, we compute the angle between the invariant manifolds at one of the transversal 
intersections. We will restrict the range of energy values to 

J G [J_,J+] = [-1.81,-1.56], (94) 

or equivalently the range of eccentricities to e € [e_,e+] = [0.48,0.67]. This is the range where 
we can validate the accuracy of our computations (see Appendix IA.4p . Below e„ = 0.48, the 
splitting size becomes comparable to the numerical error that we commit in double precision 
arithmetic. 

Remark A. 4. In this paper we concentrate on proving the existence of global instabilities in the 
Restricted three-body problem; we are not so much concerned with finding the maximal range 
of eccentricities along which the asteroid drifts. Thus we do not investigate the transversality 
of the splitting below e_. However, we are convinced that the maximal range of eccentricities 
is larger than [e_,e+], in particular that the lower bound can be pushed well below e_. We 
think that our mechanism of instability applies to this larger range of eccentricities. In fact, 
it is possible to study such exponentially small splitting using more sophisticated numerical 
methods, such as multiple-precision arithmetic, and high-order approximation of local invariant 
manifolds, see for instance [FS90l IDRR991IGS08] . 
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Figure 12: Invariant manifolds of the points p2 and p-^ for the energy level J = —1.74. 

Consider first a fixed energy level J = Jq that is close to the unperturbed situation, e.g. J = 
— 1.74. The corresponding manifolds are given in Figure [T2l In general, there are uncountably 
many intersection points. For instance, in Figure [11] we show six intersections on the symmetry 
line. However, when the perturbation is small, there is one distinguished intersection point 
located "in the middle" of the homoclinic. We call it the primary intersection point. 

Let us compute the primary intersection point zi corresponding to the "outer" splitting of 
the manifolds W''^{pz) and VF^'^fe)- For J = -1.74, the primary intersection zi corresponds 
to the first intersection of the manifolds with the = line, as we grow the manifolds from the 
fix points. Thanks to the symmetry, it is enough to look for the intersection of W'^'^{p^) with 
the = axis, because VF*'^(p2) must also intersect the axis at the same point. 

To compute the intersection point zi , we continue using a linear approximation of the local 
manifold, and propagate a fundamental domain in the local manifold by iteration. Let Vu be the 
unstable eigenvector associated to the point p^. Consider the fundamental segment between 
the points p3 + r]Vu and P{p^ + t]Vu), as in the previous section. First we look for the smallest 
natural n such that P'^[lu) intersects the = axis. Then we use a standard numerical method 
(bisection-like one-dimensional root finding) to find a point z^ in the fundamental segment 1^ 
such that 

vrp.(P"(^n))=0. 

Thus we obtain the homoclinic point zi = P'^{zu) in Figure [121 Numerically, we verify that zi 
is in the the = axis within lO^^'' tolerance. 

Finally, we vary energy J and use a continuation method to obtain the family of primary 
intersections {zi}j, using as seed the primary intersection zi\ j^_^ found above. See Figure [T3l 

Remark A. 5. For low energy levels (such as J = —1.74), corresponding to weak hyperbolicity, 
the invariant manifolds behave as if they were close to integrable, and the primary intersection 
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Figure 13: Family of primary intersection points corresponding to outer and inner splitting. 
For every energy level J, we plot the x coordinate of the intersection point zi and Z2 (the px 
coordinate is equal to zero). Notice that both families are continuous. 

corresponds to the first intersection of the manifolds with the px = axis. For high energy levels 
(such as J = —1.6), corresponding to strong hyperbolicity, the manifolds develop some folds, 
and thus the primary intersection may not correspond to the first intersection of the manifolds 
with the Px = axis. See Figure [TTJ 

In practice, we first identify the primary intersection at low energy levels, and then use a 
continuation method to obtain the primary family of intersections up to high energy levels. 

Analogously, we compute the family of primary intersections {22} j corresponding to the 
inner splitting. See Figure [13l 

Let us now compute the splitting angle between the manifolds W^'^{ps) and W^'^{p2) at 
the point zi. For illustration, we show some numerical results corresponding to the energy 
value J = —1.74. First we need the tangent vectors Wu and Ws to the manifolds at zi. See 
Figure [21 As found above, let Zu be the point in the unstable fundamental segment that maps 
to zi, i.e. P"'{zu) = z\. Consider the tangent vector Vu to the manifold ^^"'^(ps) at the point 
Zu- (Recall that at this point the linear approximation is good enough, so we can use as the 
unstable eigenvector.) Multiply Vu by the Jacobian of P at the successive iterates P^{pu), for 
i = 0,...,n — 1. This way, we obtain the tangent vector to the unstable manifold at zi. Let us 
denote this vector Wu = {wi,W2)- We normalize it to \\wu\\ = 1- 

Due to reversibility, the vector Wg tangent to the stable manifold at zi is Ws = {wi, —W2)- 
See Figure [HI Notice that we choose the tangent vectors with the appropriate orientation, i.e. 
with the same orientation as the trajectories on the manifolds. 

Thus the oriented splitting angle between Wu and Ws is 

a = 2 arctan2(— li^i, —W2), 
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Figure 14: Outer splitting of the manifolds for energy level J = —1.74. This is a magnification 
of Figure [12] at the intersection point zi. We show the vectors Wu,Ws tangent to the unstable 
and stable manifolds at zi. The splitting angle a is the angle between Wu and Wg- 



inner 


outer 


(-1.695,-1.694) 
(-1.726,-1.725) 
(-1.756,-1.755) 
(-1.781,-1.780) 
(-1.802,-1.801) 


(-1.701,-1.700) 
(-1.731,-1.730) 
(-1.760,-1.759) 
(-1.784,-1.783) 
(-1.805,-1.804) 



Table 1: Subintervals of J G [J_, J-|_] containing the zeros of inner splitting (left column) and 
outer splitting (right column). 

where arctan2 is the arctangent function of two variables, which uses the signs of the two 
arguments to determine the sign of the result. 

Finally, we let J change and, using this procedure, we are able to obtain the splitting angle 
for energy levels J G [J_, J+J. See Figure [T5l The splitting angle is nonzero for all energy values 
except for a discrete set of them. The splitting angle oscillates around zero with decreasing 
amplitude as J — t- J„ . Numerically, we find that the zeros of the splitting angle are contained 
in the intervals listed in Table [TJ 

Notice that the inner and outer splittings behave similarly. However, they become zero at 
different values of J, as seen in Table [TJ Thus, when one of the intersections becomes tangent, 
the other one is still transversal, and we can always use one of them for diffusion. 
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Figure 15: Splitting angle associated to inner and outer splitting. 



Px 






— 


-0.00002 
-0.00001 
0.00000 
0.00001 
0.00002 


-5.481541931871417 
-5.481541931790012 
-5.481541931822124 
-5.481541931967703 
-5.481541932226887 


-5.481541932226887 
-5.481541931967703 
-5.481541931822124 
-5.481541931790012 
-5.481541931871417 


0.000000000355470 
0.000000000177691 
0.000000000000000 
-0.000000000177691 
-0.000000000355470 



Table 2: Sampling of the manifolds W'^'^{p3) and W^'^{p2) at different values oi px, and their 
difference (last column). 



A. 4 Accuracy of computations 

For small eccentricities, the splitting angle a becomes very small. We need to check the validity 
of a, making sure that the size of (accumulated) numerical errors in the computation is smaller 
than the size of a. 

The smallest splitting angle in Figure [T5l corresponding to J_ = —1.81, is 

cr(J_) = -1.777970294158603 x 10~^ 

We check the validity of <t(J_) by recomputing this angle using an alternative numerical method. 
First we compute the intersection of the manifolds W^'^ips) and W'^'^{p2) with the horizontal 
axis defined by 

_ _j_ 
P- - 105 

for j E (-2,-1,1,2). 

In Table[2]we tabulate the x coordinate of VF"'^(p3) and W^'^{p2) on these axes, and their 
difference d = — gives the distance between the manifolds. We apply numerical differen- 
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Figure 16: Splitting angle (t( J) and estimate of the numerical error err( J) as a function of energy 
level J. 

tiation to the last column of this table, using central differences centered at zi with step sizes 
0.00002 and 0.00004, and obtain the values: 

_ 0(0.00001) -d(-O.OOOOl) ^ _„ „„„„i„,gi 
0.00002 

d(0.00002) - (i(-0.00002) 

d2 = — = -0.0000177735. 

0.00004 

Finally, we use Richardson extrapolation and obtain: 

d = '^^^ ~ = -0.00001776763333333333. 
Thus, using this alternative method, we obtain the splitting angle 

(t(J_) = atan(-0.00001776763333333333) = -0.00001776763333146364. 

Compare the splitting angle computed using the two methods. They differ by approximately 
10~^. This gives an estimate of the numerical error commited in our computation of the splitting 
angle. 

We repeat this test for a range of energies J € [-1.81,-1.8]. In Figure \T6\ we compare 
the splitting angle cr{J) and the estimate of the numerical error err (J). This error stays below 
10~^, and it is several orders of magnitude smaller than the splitting angle. For higher energy 
values J G [—1.8, —1.56], the splitting angle is large, so the numerical error is certainly smaller. 
Therefore we are confident that the splitting angle has been accurately computed in the range 
of eccentricities considered, [J_, J+]. 
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Figure 17: Energy J = —1.74. Resonance structure in Cartesian coordinates. The axis of 
symmetry is marked with a horizontal hne. 

B The resonance in Delaunay coordinates 

In Appendix |A] we have studied the periodic orbits and the invariant manifolds in rotating 
Cartesian coordinates {x, y,Px,Py)- Nevertheless, the study of the inner and outer maps are done 
in rotating Delaunay coordinates. Indeed since these coordinates are action-angle coordinates 
for the two body problem, it is much more convenient to use them to study the mean motion 
resonance. 

The Poincare section {y = 0} is completely different from the section {g = 0} which will 
be used from now on (see ([22]) ). In particular, the periodic orbits obtained in 

Appendix lA.ll intersect the section {y = 0} six times whereas they intersect {g = 0} seven 
times. However, we remark that the homoclinic points zi and Z2 lie on the symmetry axis both 
in Cartesian and in Delaunay variables. See Figures [T71 and [T8l 

To obtain the intersection of these periodic orbits with {g = 0} we just need to express the 
6-periodic points of the Poincare map P obtained in Appendix lA.il in Delaunay coordinates and 
then iterate them by the flow of the circular problem expressed in Delaunay coordinates until 
they hit the section {g = 0}. We do the same with the homoclinic points. In Appendix lB.il we 
explain how to compute the change of coordinates and the vector field in Delaunay coordinates. 

B.l Prom Cartesian to Delaunay and computation of dc^Hcnc 

We explain an easy way to obtain the rotating Delaunay coordinates from rotating Cartesian 
(or polar) coordinates in the circular problem. First recall that G can be computed as 

G = r {—px sin + Py cos (j)) . 
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Figure 18: Energy J = —1.74. Resonance structure in Delaunay coordinates. The symmetry 
corresponds to / = and / = vr and is marked with a vertical hne. 



The potential fiAHcnc hi Cartesian coordinates only depends on the position (x, y) of the 
asteroid, and can be easily computed. Then, one can use the equation 



J 



2L2 



G + fiAHc 



to obtain L. Knowing L and G we can obtain the eccentricity e by 



Using that r = L^{1 — ecosu), one can obtain u and from here i using Kepler's equation 
u — e sin u = i. On the other hand, from u we can obtain v using 



/ 1 + e u 

tan — = \ tan — . 

2 V 1 - e 2 

Finally, we can deduce g using that (p = v + g. 

We devote the rest of this appendix to compute dcAHdrc- The other derivatives of AHc 
can be computed analogously. Let us define 

D[ro] = D[ro]{r,v,g) = (r^ + - 2rrocos(u + g)) ^^"^ . 

Then 

AH^^rc{L,l,G,g) = -{I- ij)D[-ix]- ^lD[l- ^,]-D^. 
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Thus by the chain rule there only remains to compute dor and dev. First, let us point out that 



On the other hand, using that i = u — esinu, one has that 

sinu 



1 — e cos u 

Then, since r{L,e,i) = L^(l — e cos ii(e, using that 



we have that 
and therefore. 



cosu-e 

cost; = , (95j 

1 — e cos u 



der{L, e, i) = cos v 
Gcosv 



dGr{L,£,G) 



e 

To compute dcv, let us point out that dcv = devdce- Therefore it only remains to compute 
deV, we obtain it using formula ()95p and 



smz; 



\/l — sin u 

1 — e cos u 

Then, 



smv 
1 - e 

and therefore. 



OeV = 2 (2 + e COS V) 



smv 

Ogv = — (2 + e COS v) . 

C Numerical study of the inner and outer dynamics 
C.l Inner and outer dynamics of the circular problem 

In this appendix, we numerically compute the inner map Fq"" and the outer maps J^"*'* of 
the circular problem, given in Section [3l Recall that to compute these maps we deal with the 
extended system given by the Hamiltonian H in ()17p with cq = restricted to the energy level 
H = and thus, we have that / = — J. Then, we consider / G /-|_] = [— J+, —J~], where the 
range [— J-|_, — J_] is given in 

As seen in Section [321 the inner map has the form 

where Tj = IAtt + nTo{I) is the period of the periodic orbit obtained in Theorem [3] on the 
corresponding level of energy J, which now corresponds to an invariant hyperplane / = constant. 

Recall that we computed the periodic orbit A j as well as its period Tj in Appendix lA.li In 
particular. Figure [6] shows a plot of the function Tj — 14tt = ^Tq{I). Notice that the derivative 
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of the function To{I) is nonzero for the whole range [/_,/+]. This shows that the inner map is 
twist. Moreover, Figure E] shows that 



< fiToil) < 60fi < TT. 

Therefore, the function To{I) satisfies the properties stated in Lemma 13.31 

As a test, we have computed the same function To{I) using two different methods. First 

by computing the period of the periodic orbit, as above. Then by computing the integral 

expression (j39p using numerical integration. The difference in Toil) using both methods is of 

the order 10"^^. 

As seen in Section 13.31 the outer maps have the form 

Wf.^i.rnl' * = f'b. (97) 



t J \ t + flUJ*{I) 

For simphcity, let us only discuss the computation of uj^{I) {uj^{I) is computed analogously). 
Recall from Lemma 13.71 that the function uj^ (/) is defined as 

where, taking into account that the homoclinic orbit is symmetric with respect to the involution 

JUI) = uiil) - u;lil) = (98) 

with 

a;+ /) = Imi / — — — — ^ j— da + NTo{I) \ , (99) 

To obtain a;^(/), we compute the integrals ([99]) and (jlOOp numerically, using a standard algorithm 
from the GSL library [G"*"] . The integrals are computed within a relative error limit 10~^. 

The function dcAHdrc involved in both integrals is given explicitly in Appendix lB.il The in- 
tegral wfn(/) is evaluated on a periodic trajectory A|(a) of the reduced circular problem (namely, 
with reparameterized time, see (j3ip ) with initial condition ^4, a fixed point of the Poincare map 
Vq found in Appendix|Bl The integral w^j_(/) is evaluated on a homoclinic trajectory 7/(cr) of the 
reduced circular problem with initial condition Z2, the primary homoclinic point corresponding 
to the inner splitting found in Appendix IA.31 

Next we make a couple of important remarks about the numerical computation of the integral 
u!^_^{I). The key point is that the homoclinic orbit 7^ was already computed in Appendix IA.3I 
with high accuracy, and we can exploit this information here. Recall that the primary homoclinic 
point Z2 was obtained as the n-th iterate of a point Zu in the local fundamental segment under 
the Poincare map: 

Z2 = {viriz^). (101) 

Moreover, recall that the point Zu was chosen to be suitably close to the fixed point for each 
energy level J. See Remark lA. 31 

Notice that the integral w^j_(I) is defined by a limit as — )• 00, i.e. as the homoclinic orbit 
7}(cj) asymptotically approaches the periodic orbit Aj(cj) in forward time (see equation ([10]) ). 
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Figure 19: Functions and involved in the definition of the outer map (j97p of the 

circular problem as a function of the Jacobi constant J (recall that in the circular problem 
I=-J). 

Numerically, of course, we should stop integrating at an upper endpoint N large enough such 
that the integral converges. In practice, we choose the upper endpoint N = N{I) to be the 
number of iterates n = n{I) in (|10ip . This means that we evaluate the integral along the 
homoclinic trajectory 7|(o") until it reaches the point z^, which is suitably close to the periodic 
orbit. 

Notice also that integrating the homoclinic trajectory 7|(o") forwards in the reduced system 
means integrating it backwards along the unstable manifold in the original system. This is 
numerically unstable, since numerical errors grow exponentially. In practice, we rewrite the 
integral ()99p using the change of variables a = a — 14A^7r so that the homoclinic trajectory is 
integrated forwards along the unstable manifold, starting from the point 

The computed values of the functions and are shown in Figure [T9l Note that 

they are plotted as a function of the Jacobi constant J instead of as a function of /, so that they 
can be compared with Figure [6l where we have plotted ^To{I) = Tj — lAn as a function of J. 

To test the computation of the function uj^^, we directly verify the definition of the outer map 
in 13.41 Let Z2 = {Li^,lii,Gh,Q) be the primary homoclinic point, and let ps = {Lp,£p,Gp,0) be 
the periodic point. Given a point {Lh,ih,Gh,0, 1,t) in the extended circular problem, we check 
that it is forward asymptotic (in the reparametrized time) to the point (Lp, ip, Gp, 0, /, t+a;^j_(I)), 
where t G T is arbitrary. Thus we check that the distance 

dist+(s) = \<^o{s, {Lh, ih, Gh,0, 1, t)} - <^o{s, {Lp, ip, Gp, 0, /, t + a;^(/))}| 

with exponential decay. 

The result of the test is shown in Figure for values of the energy J G [J_, J+] (recall that 
J = —I). Notice that the vertical axis is in logarithmic scale. Let s = 14A^7r. We plot the 
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Figure 20: Exponential decay of the function dist"*" as a function of N (multiples of the period) 
for different energy levels. The energy levels J G [ J_ , J+] are color-coded. 



distance dist"*" as a function of (multiples of the period). The test shows exponential decay 
of the distance function for all energy values, i.e. straight lines in the plot. 

Recall that the periodic orbits Aj^(s) become more hyperbolic as the energy / decreases. 
Thus, the rate of exponential convergence between the homoclinic and the periodic trajectory 
also increases, i.e. the straight lines have increasing slope in the plot. As explained above, the 
length of integration = A^(/) along the homoclinic orbit is suitably chosen for each energy 
level. For /—)■/_, there is exponential decay up to time s = 40 • (IAtt) ~ 1760. 



C.2 Inner and outer dynamics of the elliptic problem 

In this appendix, we numerically compute the first orders in cq of the inner map J^" and the 
outer maps 7-"°"*'* of the elliptic problem, given in Section [H In order to compare the inner 
and outer dynamics of the elliptic problem through Lemma |5.2|, only some specific terms in the 
expansions of the inner and outer maps are necessary. Namely, we only need to compute the 
term Ai in the expansion of the inner map ([62]), and the terms B* in the expansion of the outer 
maps ([75]) . 

Recall from Section [4.41 that Ai can be split as 

Ai{I,t)=At{I)e'' + A^{I)e-''. 

Since Af and A^ are complex conjugate, it is only necessary to compute one of them. Let us 
compute the positive harmonic, 

Jo -1 + l^ldGAHcircliicr) 
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Figure 21: Function Af{I) (real and imaginary parts) involved in the definition of the inner 
map ()62p of the elliptic problem as a function of the energy of the system in rotating coordinates 
H. Recall that H = -I. 



Notice that the denominator is the same one used in the previous section for the inner and outer 
dynamics of the circular problem. Next we give the numerator iAH^^^ explicitly. Let 



AH^,,{L,£,G,g,t) 



1-^ 



r{L,£,G) 
r{L,e,G) 



,v{L,£,G),g,t] , 



1 — fj. \ 1 — 

where Bi is the function defined in Lemma 14.41 Then it is straightforward to see that 



AHl^{l,L,g,G) 



Bt 



-Bi 



r{L,i,G) 
A* 

r(L,AG) 



,viL,i,G),g 



viL,i,G),g 



where 



Bt{r,v,g) 



1 — r cos(f + g) — i2r sm{v + g) 
2A^{r,v + g) 



(103) 



The computed value of the function Af is shown in Figure [211 We plot it as a function of 
the energy of the elliptic problem in rotating coordinates H in ()16p . Recall that since we are 
working in the energy level H = of the extended Hamiltonian H in (|17p . we have that / = —H. 

For the outer map, we compute the functions B*{I). Similarly to Ai, it is only necessary 
to compute the positive harmonics B*''^. Recall from Lemma 14.91 that the positive harmonics 
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Figure 22: Functions B^'~^ and B^'~^ (real and imaginary parts) involved in the definition of the 
outer map ([73]) of the elliptic problem. 



B^'+{I) and B^'+{I) are defined as 

i?f'+(/) = i?it(/) + i?rU)e 
B'''+{I) = B^-'{I) + Btt{I)e 



(104) 



where and u)\^ were obtained in Appendix IC.ll To obtain B*^'^ and -Bin"*", we compute 
the integrals (f75]) - ([77|) numerically, using the same techniques as in the previous Appendix lC.il 
In particular, the integrands of the Melnikov integrals (|75p and (|76|) . by construction, decay 
exponentially as T — >• ±00 and we take the same approximate limits of integration ibl47rA'^ 
where N = N{I) is the constant considered in Appendix lC.il 

The computed values of the functions B^'^{I) and B^'^{I) are shown in Figure [ 



C.3 Comparison of the inner and outer dynamics of the elliptic problem 

Finally, we verify the non-degeneracy condition 

B*^^ (X) / for X e (105) 

stated in Lemma 15.21 which implies the existence of a transition chain of tori. Since S*'"*" and 
B*'^ are complex-conjugate, it is only necessary to compute one of them. Let us compute the 
positive harmonic, 

B*'^ (I) = B*'^ (I) - -^^^^At {1) . 

All the functions involved in the expression above are known: To and uj* are obtained in Ap- 
pendix [CTT] and and are obtained in Appendix I C. 2 1 
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Figure 23: Functions _B^'+ and (real and imaginary parts) 



The computed values of the functions and i?^'"*" are shown in Figure [23l Therefore, 
we see that the functions i?*'"*" are not identically zero. This justifies the statement ()89p in 
Lemma I5.2i 

Remark C.l. Figure [23] also shows that and B^'~^ are almost identical, which is surprising 
for the authors. However, this fact is not relevant for the argument in Lemma [521 we only need 
that these functions do not vanish identically. 



D Conjectures on the speed of diffusion 

D.l Speed of diffusion for a priori unstable systems and Positive measure 

Consider the following nearly integrable Hamiltonian system proposed by Arnold |Arn64j : 

Hsip,q,I,(l),t) = ^p^ + cosq-l + h^ + eHi{p,q,I,4>,t), wherep,/GM, (p,q,teT (106) 

for an analytic perturbation eHi. This system is called a priori unstable. Proving Arnold 
diffusion for this system consists in showing that, for all small e > and a generic eHi, there 
exists orbits with 

\m-m\>o{i), 

where 0(1) is independent of e. 

There has been fascinating progress in this problem achieved by several groups (see jBerOSl 
l(]Y04[ iDdlLSMl IDH09I ITren4j l. Treschev |Tre04j not only proved existence of Arnold diffusion, 
but also gave an optimal estimate on speed, namely, he constructed orbits 



|/(i) - 1(0)1 >c. 
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for some c > 0. One can see that this estimate is optimal, i.e. \I(t) — /(0)| < C 



e 



t for some 



Ine 



C>c. 



Heuristicahy the mechanism of diffusion is the fohowing. For small e > the Hamiltonian 
iJg has a 3-dimensional normally hyperbolic invariant cylinder Ag close to Aq = {p = q = 0}. A 
hypothetical diffusing orbit starts close to and makes a homoclinic excursion. Each homoclinic 
excursion takes approximately 0{\ lne|)-time. Increment of I{t) after such an excursion is ©(e)!! 
If one can arrange that all excursions lead to increments of I{t) of the same sign, the result 
follows. 

It seems natural that orbits will be trapped inside the resonance p = for polynomially long 
time. Using this heuristic description one can conjecture that increments I{t) can behave as a 
random walk for positive conditional measure for polynomially large time. 

Positive measure conjecture H Consider the Hamiltonian with a generic perturbation 
eHi. Pick an e-ball of initial conditions, whose center projects into {p,q) = 0, and denote 
the Lebesgue probability measure supported on it by Lebg. Then, for some constants c, C > 
independent of e, the set of initial conditions satisfying 



Since a typical excursion takes 0(1 lne|)-time and each increment is 0{e), we essentially 
conjecture that after excursions with uniformly positive probability there will be drift 

of order 0{e)0{e-^) = 0{l). 

D.2 Structure of the restricted planar elliptic three-body problem 

In this appendix we relate a priori unstable systems and the elliptic restricted planar three-body 
problem. Recall that we managed to write the Hamiltonian of the latter problem in the form 

Hc\\{L, e, G, g, t) =HciTc{L, i, G, g, /u) + fieoAHcn{L, i, G, g, t, fi, cq) 



We have that 

• Hq is an integrable Hamiltonian. 

• Hcirc is non-integrable and for Hdrc in a certain interval [J_, J+] there is a family of 
hyperbolic periodic orbits {pj} whose invariant manifolds intersect transversally along at 
least one homoclinic. 

• -ffeii is a 0(/ieo)-perturbation of Hdrc such that a Melnikov integral evaluated along a 
transverse homoclinic being non-degenerate in two different ways: dependence on time is 
non-trivial and relation between inner and outer integrals is non-degenerate (see ([9])). 

^This is only an heuristic description as dynamics inside of the cyhnder should come into play. Near so-called 
double resonance dynamics is different from the one near single resonances 

® This conjecture is in a sense answers a question posed by Katok about measure of diffusing orbits 



|/(r)-/(o)|>i 



for some 



< T < C 



Ine 



is denoted Diff and has measure Leb ^ ( Diff) > c. 
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Having all these non-degeneracy conditions we prove the existence of diffusing orbits. 

It is not difficult to prove, using averaging techniques, that for ^ > small, there is a family 
of saddle periodic orbit {77} on some interval [J_, J+]. It seems, however, to be a non-trivial 
problem to establish the splitting of its separatrices. Due to reversibility ([HI]) there are at least 
four homoclinic intersections (two for upper separatrices and two for lower ones). Having these 
two conditions it is natural to expect that at least one of the four associated Melnikov integrals 
is non-degenerate. Qualitative analysis shows that it should be possible to have a homoclinic 
excursion 0(//eo)-close to the invariant cylinder. Such an excursion takes 0{\ ln(/zeo)|)-time. If 
the excursion is selected properly, then the result of the excursion is that the increment of the 

eccentricity is O{fieo). This makes us believe that the instability time obeys T ~ ttt^ 

stated in ([6]). 

Let us point out that we believe that our diffusion mechanism survives even for non- 
infinitesimal eo's, e.g. realistic eo = 0.048. To justify this, we review the above structure. 

Notice that we use a 3-dimensional normally hyperbolic invariant cylinder and the intersec- 
tion of its invariant manifolds to diffuse. The cylinder arises from the family of periodic orbits 
{7j} Je[J-,J+] of the circular problem, which persist under the elliptic time-periodic perturbation 
fieoAHcii{L, £,G, g,t, HjCq) (see (fTUj) and the derivation in the corresponding section). As the 
analysis carried out in Section [4.21 shows, in the neighborhood of this family {7j}j6[j_,j^], the 
perturbation fieoAHgii{L,£,G, g — t,t,fj.,eo) can be averaged out to ©(^efj). Thus, invariant 
cylinders could persist even for not very small eo's. However, estimating remainders analytically 
after several steps of averaging is nearly impossible. Numerically though it might be feasible. 

Once the existence of an invariant cylinder is established, we need to justify the existence 
of transverse intersections of its manifolds. As before, analytically it is an insurmountable task, 
but numerically it seems to be an achievable goal. 

If these two steps are done, then one could try to compute numerically inner and outer 
maps and show that they do not have common invariant curves. This is again a difficult, but 
numerically realistic task (see |DMR08] for the computation of the outer map in another problem 
in Celestial Mechanics). 

On the other hand, the above asymptotics probably does not hold in the neighborhood of 
circular motions of the massless body, which might be much more stable than more eccentric 
motions. Yet many other factors might influence the local stability or instability of various 
objects (see section [1.2.3p . 



D.3 The Mather accelerating problem and its speed of diffusion 

The structure we use to build diffusion is similar to the Mather acceleration problem. Let us 
recall this problem and state an interesting result of Piftankin jPif06j on speed of diffusion. 
Consider a Hamiltonian system 

Hiq,p,t)=K{q,p) + V{q,t), q e T\ p G M^, t € T, 

where K{q,p) = ^{A~^{q)p,p) — kinetic energy corresponding to a riemannian metric K{q,p) = 
^{A~^{q)p,p), p = A{q)q, q G TqT^ and V{q,t) is a time-periodic potential energy. Since the 
system is not autonomous energy is not conserved. 

HI Suppose the geodesic flow associated to K has a hyperbolic periodic orbit T and transversal 
intersection of its invariant manifolds, which contains a homoclinic orbit 7(t), t G M. 
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H2 The Melnikov integral is not constant. More exactly, define a function 

/T f-T+t' 
V{-f{t),t)dt- / V{-f{t),t)dt. 

The limit turns out to exist and is independent of a choice of t^, t^. This function is 
assumed to be non-constant. 

Mather and his followers |Mat96[ IBT991 IDdlLSOni \GT08[ IKalOSI iPifOH] proved existence of 
an orbit {qT{t),PTit)), t G M of unbounded energy, de la Llave |dlL04) . Piftankin |Pif06j . and 
Gelfreich-Turaev |GT08] proved that such an orbit can be chosen to have linear growth of energy 

H{qr{t),PT(.t)) >At + B for all t > 

for some ^ > and G M. 

Notice that for large energies H ~ the conformal change of coordinates 

p = -, H = e^'^H, t = ei 

£ 

leads to the new Hamiltonian 

= i^(g,p) + eV(g,ef). 

It was shown in |dlL04l IPifOGl IGTOSj that there are orbits diffusing linearly in the size of the 
perturbation. In order to see these orbits, notice that K{q,p) has a horseshoe. Then, H can 
be considered as a time-periodic perturbation over such a horseshoe. It is shown by different 
methods in |dlL041 IPif06l IGTOSj that for a generic time-periodic perturbation of a horseshoe 
there are linearly diffusing orbits. 

D.4 Modified positive measure conjecture 

For systems with the properties discussed above we can modify the positive measure conjecture 
as follows: 

Positive measure conjecture for Mather type systems Consider the Hamiltonian 

H^^e{L, £, G, 5, t) = HI (L, G) + fiAHo{L, £, G, g, fi) + fieoAHi{L, I, G, g, t, fi, cq) 
such that 

• for some interval [J_, J+J the Hamiltonian Hq + [iAHq has a family of saddle periodic 
orbits {pj}jG[j_,j+], 

• for each J G [J_J+] there is at least one transverse intersection of its invariant manifolds, 

• A Melnikov integral evaluated along a transverse homoclinic and inner dynamics are non- 
degenerate: the dependence of the Melnikov integral on time is non-trivial and the relation 
between inner and outer maps is non- degenerate (see 
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Pick a fieQ-ball of initial conditions B^^o whose action components are centered at a resonance 
between i and g. Denote the Lebesgue probability measure supported on the ball B^^^ by Leb. 
Then for some constants c, C > independent of /U and cq, the set of initial conditions satisfying 

|G(r)-G(0)| > 1 for some < T < C^^^ 
is denoted Diff and has measure Leb ( Diff) > c. 



Here is an important difference between the system -ff^,eo priori unstable one H^, 

given by ()106p : the Hamiltonian Hq + fiAH^ already has "chaos" and a family of horseshoes 
on each energy surface with J G [J_ J+], while Hq = — eHi is integrable. As we pointed 
out above, for a generic time-periodic perturbation over a horseshoe there are orbits diffusing 
linearly fast |dlL04l |Pif061 lUTOSj . Yet we are interested in a set of conditional positive measure. 

In order to see the time of diffusion on an heuristic level, notice that Hq + ^iAHq has a 
family of saddle periodic orbits {pj}j£[j_.j^] whose exponents are ~ -^//I. Thus, one homoclinic 
excursion passing //eo-close to separatrices takes | ln/ieo|/y^-time. Each excursion might lead 
to increment of G of size ~ /ieo- Conjecturing that random walk approximation holds true to 
have C'(l)-changes in G, we need 0{fi^'^eQ^) excursions. 
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